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1.1  OVERVIEW  OF  PROGRAM 
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1_  Problem  Definition 

The  problem  addressed  in  this  program  is  that  of  processing  information 
produced  at  the  output  of  an  infrared  detector  array  to  determine  the  pres- 
ence or  absence  of  intruders  in  the  field  of  view  of  each  element  of  that 
array.  The  change  produced  by  each  detector  in  the  array  is  stored  in  a 
charge-coupled  device  (CCD)  and  read  out  serially  at  periodic  intervals. 
The  physical  and  electrical  characteristics  of  the  detector  array  and  CCD 
are  assumed  to  be  known. 

It  is  convenient  to  divide  the  problem  into  three  major  area.  These 

are: 


a)  What  is  the  response  of  the  detectors  to  intruders  having 
specified  size,  speed  and  thermal  properties? 

b)  What  signal  processing  techniques  are  most  suitable  for  dis- 
tinguishing between  responses  created  by  intruders  of  interest 
and  response  created  by  irrelevant  objects? 

c)  What  performance  can  be  achieved  with  the  candidate  system 
that  appears  to  most  nearly  meet  the  requirements? 

These  three  aspects  of  the  problem  are  discussed  in  detail  in  the  following 
sections. 
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1_.  1_.£  Objectives  of  the  Investigation 

The  objectives  of  this  investigation  follow  immediately  from  the  prob- 
lem definition  stated  above.  However,  a more  precise  statement  of  these  ob- 
jectives is  desirable  as  a means  of  both  summarizing  the  scope  of  the  inves- 
tigation and  indicating  the  sequence  in  which  the  problems  were  addressed. 

Thus,  the  objectives  may  be  described  as  follows: 

1)  Develop  reasonable  models  for  the  specified  classes  of  in-  ► 

truders  and  backgrounds. 

2)  Determine  the  detector  response  to  each  of  the  intruder  and 
background  models. 

3)  Carry  out  a theoretical  analysis  of  the  basic  detection  prob- 
lem. J 

4)  Conceive  various  practical  implementations  of  the  detection 
scheme  and  evaluate  their  relative  merits. 

5)  Select  a detection  scheme  and  identify  the  optimum  parameters 
for  that  scheme. 

6)  Evaluate  the  performance  of  the  selected  detection  method  with 
respect  to  the  specified  intruder  models. 

7)  Prepare  a preliminary  design  and  rough  cost  estimate  for  a 
breadboard  realization  of  the  proposed  detection  method. 

• I 

1^.1_.3^  Brief  Summary  of  Results 

Although  the  main  body  of  this  report  is  a detailed  presentation  of  the 
results  of  this  investigation,  a brief  summary  of  the  major  achievements 
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provides  a useful  introduction  to  this  material.  All  of  the  objectives  out- 
lined above  have  been  "'•hieved. 

1)  Intruder  models  have  been  developed  for  alt  of  the  specified  classes. 
These  models  are  simple  enough  to  permit  computer  evaluation,  but  also 
represent  accurately  the  size,  speed,  temperature  and  thermal  emissivity 
of  each  class  of  intruder.  Models  of  the  background  have  also  been 
developed  for  a wide  range  of  circumstances.  In  addition,  the  effects  of 
solar  illumination  and  atmospheric  precipitation  have  been  modeled  on  a 
simpl if ied  basis. 

2)  A computer  program  has  been  developed  that  computes  the  detector  output, 
as  a function  of  time,  as  any  intruder  model  enters  the  field  of  view  at 
any  angle  and  with  any  speed. 

3)  The  theoretical  optimum  detector  has  been  analyzed  and  evaluated  for 
cases  that  permit  comparison  with  the  practical  implementation. 

A)  A proposed  practical  implementation  has  been  selected  and  the  optimum 
parameters  determined  by  "worst  case"  analysis. 

5)  The  proposed  detection  method  has  been  evaluated  with  respect  to  the 
background  and  intruder  models  by  calculating  the  probabilities  of  false 
alarm  and  the  probabilities  of  detection  for  each  of  the  models.  Two  im- 
portant conclusions  from  this  study  are: 

(a)  The  probabilities  of  false  alarm  due  to  system  and  background 
noise  are  negligibly  small  when  the  decision  thresholds  are 
adjusted  to  reject  non-threat  classes  of  intruders. 
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(b)  The  probabilities  of  detection  for  other  classes  of  intruders 
are  essentially  unity  under  the  same  threshold  conditions. 
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The  major  deficiency  of  the  proposed  system  appears  to  be  the  possibility  of 
false  alrms  due  to  scattered  clouds  or  sun  glint. 


6)  The  preliminary  design  of  a mi croprocessor  that  can  accomplish  the  neces- 
sary data  processing  for  the  detection  algorithm  has  been  prepared.  The 
estimated  cost  of  the  electronics  components  only  for  a breadboard  model 
on  a "make-one"  basis  is  about  $1200  and  the  dc  power  requirements  are 
about  15  amperes  at  5 volts. 


7)  It  is  recommended  that  the  development  of  a working  model  be  initiated 
and  some  steps  to  be  taken  in  this  direction  are  outlined. 


2.1  SIGNAL  CHARACTERISTICS 

1_.1_  Detector  Response 

The  intent  of  this  section  and  the  two  following  sections  is  to  state 
the  principal  equations  that  have  been  used  to  calculate  the  signals  that 
are  generated  by  the  detector  array.  These  equations  have  been  obtained 
from  the  literature  and  from  information  supplied  by  ETSD  personnel  and  no 
attempt  has  been  made  to  attribute  them  to  specific  sources. 

The  detector  elements  are  Schottky  interval  emission  photodiodes  havi>>g 
an  electron  yield  of 


(v)  = 


n v v 


(•electrons. 
v photon 


(1-1) 


where 


/ 


V 

W 


h = 6.6256  x 10  ^ (Ws^)  (Planck's  constant) 


v = photon  frequency  (s  ) 


* = barrier  height  (Ws) 


and  is  a factor  determined  by  the  geometrical,  optical  and  transport  pro- 
perties of  the  photodiode.  It  is  convenient  to  express  the  electron  yield 
in  terms  of  wave  number,  v,  since  the  eventual  integration  is  carried  out  in 
this  variable.  Thus,  if 


0 = * 


(1-2) 


he (C.j  (v-\Jq) ) 


(1-3) 


where 


"0  = ^ms/hc 


c = 3 x 10,u  cm/s 


For  a black  body,  at  absolute  temperature  T,  the  spectral  radiant  pho- 
ton emittance  as  function  of  wave  number  is 


♦ (v,T)  = 


exp(vhc/kT)-1 


(1-4) 


where 
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k = 1.3 8 x 1(f23  Ws/°K  (Boltzmann's  constant) 


2 -1 

in  photons  per  second,  per  steradian,  per  cm  , per  cm  wavenumber.  If  the 
black  body  is  at  a distance  R from  the  detector,  and  if  it  is  imaged  onto  a 
detector  cell  by  lense  having  a diameter  D,  then  the  solid  angle  subtended 
by  tha  detector  is 


fl  = O^)2  (steradian) 


(1-5) 


Furthermore,  if  the  target  has  an  area  of  Ag(cni  ) and  its  image  completely 
covers  a detector  cell  having  an  area  of  Ag(cm2),  then  the  number  of  photons 
arriving  at  the  detector  cel l is 


wA  <3 
eJ 

2F2 


4>  ( v,T) 


in  photons  per  second  per  cm  1 wavenumber  where  J is  the  transmittance  of 
the  optics  ( £ 1)  and  F is  the  f-number  of  the  optics. 

The  number  of  electrons  produced  in  a stare  time  t when  the  detector 
element  is  filled  by  a black  body  at  temperature  T can  now  be  expressed  as 


*A3t  2 
e s 


BB 


2F 


j ♦(v,T)n(v)dv  (electrons) 


(1-6) 


where  and  are  determined  by  the  detector  cutoff  and  the  optical  system 
cutoff. 

The  electrons  produced  by  the  detector  are  converted  to  voltage,  stored 
in  the  CCD,  and  subsequently  read  out  as  voltages.  However,  under  the  as- 
sumption that  relationship  between  electrons  and  voltage  is  a linear  one,  it 
does  not  really  matter  what  the  constant  of  proportionality  is  since  all  of 
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the  subsequent  computations  are  performed  in  terms  of  ratios. 

The  number  of  electrons  given  by  (1-6)  is  in  reality  the  mean  value  of 
the  number  of  electrons  collected  in  one  stare  time.  The  actual  number  of 
electrons  in  each  stare  time  is  a random  variable  having  a Poisson  distribu- 
tion. One  of  the  characteristics  of  the  Poisson  distribution  is  that  the 
mean  and  variance  are  identical  Cl,  p.  1 453 . However,  the  Poisson  distribu- 
tion, being  discrete,  is  difficult  to  handle  analytically  in  making  studies 
of  probability  of  detection  and  probability  of  false  alarm.  Fortunately, 
the  number  of  electrons  is  so  large  (greater  than  10"*)  that  the  Gaussian  ap- 
proximation to  the  Poisson  distribution  is  extremely  good.  Thus,  the  elec- 
trons produced  in  each  stare  time  will  be  assumed  to  be  a Gaussian  random 
variable  with  a mean  and  variance  given  by  (1-6),  or  by  modifications  of 
(1-6)  to  be  discussed  subsequently. 

£.1^.£  Background  Signal 

In  the  absence  of  an  intruder  each  cell  of  the  array  is  viewing  only 
background.  In  order  to  calculate  the  signal  produced  by  this  background  it 
is  necessary  to  make  two  modifications  to  (1-6).  In  the  first  place,  the 
background  is  not  a black  body.  For  all  of  the  calculations  made  here,  it 
will  be  treated  as  a gray  body  with  an  emissivity  of  e0  (0  < e0  £ 1.0)  that 
is  not  a function  of  wavelength.  Thus,  it  is  only  necessary  to  multiply 
(1-6)  by  £g  to  account  for  this. 

Secondly,  it  is  necessary  to  account  for  attenuc  ion  of  the  background 
emitted  signal  by  the  atmosphere.  This  attenuation  is  described  by  the  at- 
mospheric transmittance  defined  as 

-a.R 

t(R)  * e 1 (1-7) 
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where  a.,  is  the  attenuation  coefficient  in  nepers/meter  and  R is  the  dis- 
tance to  the  background  in  meters.  It  is  customary  to  express  attenuation 
coefficients  in  dB/km  but,  because  of  the  short  distances  involved  here,  a 
value  in  dB/m  is  used.  Thus,  let 


10a. 

= 4.343a.,  (dB/m) 


so  that 


: (R) 


. 2302aR 


(1-8) 


The  original  computations  were  performed  with  tabulated  data  for  t(R) 
that  were  wavelength  dependent.  However,  because  of  the  limited  number  of 
conditions  for  which  data  was  available,  and  because  the  effect  of  atmos- 
pheric attentuation  was  minimal,  it  was  decided  to  employ  (1-8)  instead. 
Values  of  a rang'ng  from  0 to  .07  dB/m  appear  to  cover  all  conceivable  at- 
mospheric conditions. 

It  is  further  assumed  that  the  atmosphere  is  emitting  radiation  as  a 
black  body  at  the  same  temperature  as  the  background.  Thus,  the  numbe>-  of 
electrons  resulting  from  background  emission  at  temperature  T0  is 


nAOt  52 

N0E  = 2 - f "Ce  B<t>  (v,Tg)x  (R)  + *(v,TB)C1-T(R)]>n<v)dv  <1-9) 

2F  *1 

The  result  given  by  (1-9)  is  applicable  for  both  day  and  night  opera- 
tion but  during  the  daytime  there  is  an  additional  signal  resulting  from  the 
reflected  signal  from  background.  There  are  many  possibilities  for  this  but 
the  only  one  considered  here  assumes  that  the  reflection  is  due  entirely  to 
direct  sunlight.  In  this  case  the  electrons  resulting  from  the  reflected 
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signal  only  becomes 


"A3t  v2 
e s f 


V (1‘tB)(^)cosesun  ♦(v,Tsun).(R)n(3)dS 


(1-10) 


where 


Asun  = area  °*  the  sun 


Rsun  = distance  to  the  sun 


0sun  = an9*e  incidence  on  the  background  surface 


Tsun  = 6000°K'  temperature  of  the  sun. 

Because  of  the  dependence  on  sun  angle,  all  calculations  were  made  with  a 
compromise  value  of  esun  = 30°. 

On  the  basis  of  the  above  discussion,  the  background  signal  may  be  ex- 
pressed as: 


Night:  N0  = Nbe 


Day:  NB  = NBE  + NBR 


(1-11) 


(1-12) 


Some  of  the  parameters  that  enter  into  the  computation  of  these  values,  and 
were  held  fixed  throught  the  computations,  are: 

Ae  = 5.16  x 10'9  cm2  (1  mil  x 8 mil) 
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*3  = 

F = 1.0 

C1  = 0.1  (ev)"1 

Unless  otherwise  noted,  the  wave  number  limits  were  taken  to  be 
v,j  = \>q  = 2380  (X  = 4.2  pm) 


v2  = 2945  (X  = 3.4  ym) 


2^2*!  Target  Signal 

When  an  intruder  (target)  enters  the  field  of  view  of  any  cell  there 
will  be  a change  in  the  photon  arrival  rate.  This  change  will  depend  upon 
the  emissivity  of  the  target  relative  to  the  background,  the  temperature  of 
the  target  relative  to  the  background,  and  the  fraction  of  the  field  of  view 
of  that  cell  that  is  occupied  by  the  target. 

The  number  of  electrons  that  would  be  produced  by  emission  from  a tar- 
get that  fills  the  entire  field  of  view  would  be,  by  analogy  to  (1-9), 

*A9t  ”2 

Nje  * — f <eT^(;,TT)i(R)  ♦ ♦(v,TT)C1-T(R)D>n(v)dv  (1-13) 

2F  31 

* . 

Similarly,  in  the  daytime  there  will  be  an  additional  number  of  electron  due 
to  the  reflected  signal  from  the  sun.  This  is,  by  analogy  to  (1-10), 
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(1-14) 


*A‘Jt  v2 

A 

e s c 

(1  w sun 

2F2  / 

u ETM  2 

R* 

V1 

sun 

In  both  of  the  above  expressions  eT  is  the  emissivity  of  the  target  and  TT 
is  the  temperature  of  the  target. 

When  the  target  does  not  fill  the  field  of  view  of  the  cell  in  ques- 
tion, the  above  numbers  must  be  modified.  Specifically,  if  the  target  occu- 
pies a fraction  6 (0  < i <_  1)  of  the  field  of  view,  the  target  signal  may  be 
expressed  as 


Night:  N_  = «NTC  + (1-«)ND1_ 

1 1 bt 


(1-15) 


Day:  Ny  = 6(NJE  + NJR)  + (1-6) (NQE  + N0R) 


(1-16) 


2^.4^  Signal  -to-Noise  Ratio 

The  most  important  parameter  in  determining  the  detectability  of  any 
target  is  the  signal -to-noi se  ratio.  This  is  defined  as  the  ratio  of  the 
square  of  the  change  in  mean  value  due  to  the  target  to  the  variance  of  the 
signal  when  there  is  background  only.  Since  the  mean  value  and  variance  are 
identical,  this  signal -to-noise  ratio  may  be  expressed  as 


SNR  = 


(nt  - V 


(1-17) 


Under  conditions  of  both  night  and  day,  this  becomes 
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Night:  SNR  = 


(1-18) 


6 (Nte  " nbe)2 


N 


BE 


Day:  SNR  = 


6 C(NTE+NTR)  - (NBE+NBR):|2 


nbe  + nbr 


(1-19) 


2.2  SYSTEM  DESIGN  CONSIDERATIONS 

2. 2. 1_  Signal  Processing  Requirements 

The  basic  objective  of  the  signal  processing  is  to  determine  when  the 
mean  value  of  the  signal  produced  by  each  detector  cell  changes  from  its 
normal  value  by  an  amount  sufficient  to  indicate  the  presence  of  a target  in 
the  field  of  view  of  that  cell.  Obviously,  in  order  to  accomplish  this  ob- 
jective it  is  necessary  to  establish  what  the  normal  value  is.  Unfortunate-  J 

> 

ly,  the  normal  value  may  be  different  for  each  cell  and  it  will  probably  be 
a slowly  changing  function  of  time.  Because  of  the  large  number  of  cells  in 

ft 

each  array,  this  operation  alone  represents  a substantial  computational 
load.  However,  there  are  other  operations  that  must  be  performed  regardless 
of  what  technique  is  selected  to  make  a decision.  Hence,  the  purpose  of 
this  section  is  to  list  some  of  these  signal  processing  requirements  as  a 
preliminary  to  selecting  a proposed  system. 

A/D  Conversion:  Because  of  the  large  number  of  operations  that  need  to  be 
performed  it  appears  that  the  only  feasible  approach  is  to  accomplish  them 
digitally.  Thus,  the  first  step  in  any  system  is  to  convert  the  analog  sig- 

m w 

nals  that  come  from  the  CCD  into  digital  form.  The  major  considerations  in 
performing  this  conversion  are  the  number  of  quantizing  levels  used  and  the 
amplitude  separation  represented  by  these  levels.  There  is  little  point  in 


I 
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having  each  quantum  smaller  than  the  rms  noise  associated  with  the  signal 
from  each  cell.  It  is  anticipated  that  the  rms  noise  will  be  on  the  order 
of  1 millivolt. 

The  typical  mean  value  of  the  voltage  from  each  cell  will  range  from 
100  to  400  millivolts,  and  can  be  limited  to  this  range,  or  a smaller  range, 
by  controlling  stare  time.  Hence,  it  appears  that  256  amplitude  levels, 
separated  by  1 millivolt,  should  provide  an  adequate  range.  Furthermore, 
256  levels  requires  only  8 bits  out  of  the  A/D  converter  and  this  is  a con- 
venient number  for  most  microprocessors  to  handle. 

Since  the  CCD  is  read  out  serially,  a single  A/D  converter  is  all  that 
is  required  for  each  array.  Each  conversion  should  be  accomplished  in  about 
4 microseconds. 

Stare-Time  Cont rol ; The  photon  count  on  each  cel l of  the  array  may  change  by 
factors  of  50  to  100  between  day  and  night  operation  or  because  of  changes 
in  background  conditions.  It  appears  to  be  necessary,  therefore,  to  be  able 
to  change  the  stare  time  in  order  to  avoid  CCD  saturation  at  one  extreme  or 
signals  at  the  noise  level  at  the  other  extreme.  Because  the  ambient  il- 
lumination affects  all  cells  in  the  array,  a reasonable  approach  to  develop- 
ing an  appropriate  control  signal  is  to  determine  an  average  signal  over  all 
cells,  or  at  least  over  a subset  of  cells  that  spans  the  overall  field  of 
view.  Hence,  another  operation  that  must  be  performed  in  any  system  is  that 
of  averaging  signals  across  the  array  at  any  one  stare  time. 

Reference  Signal  Construction:  Any  decision  operation  must  be  based  on  com- 
paring the  most  recent  output  from  the  array  with  a reference  signal  of  some 
sort.  There  are  a variety  of  possibilities  for  the  choice  of  a reference 
signal,  including: 
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a)  Previous  array  output 


b)  A time  average  of  previous  array  outputs 

c)  Adjacent  celt  outputs  at  the  same  time 

d)  A combination  of  (a)  and  (c). 

Regardless  of  which  method  is  selected,  however,  some  operations  must  be 
performed  on  the  signals  from  the  array  in  order  to  construct  the  desired 
reference  signal . 

Comparison  Operation:  Several  different  methods  of  making  the  comparison 
between  the  reference  signal  and  the  current  array  output  are  also  possible. 
The  most  common  operations  are  subtraction  and  correlation. 

Decision  Operation:  After  the  comparison  has  been  performed,  it  is  necessary 
to  make  a decision  as  to  the  significance  of  the  comparison  result. 
Although  the  ultimate  decision  is  binary,  (i.e.,  either  there  is  an  intruder 
or  not)  that  decision  may  be  arrived  at  by  making  several  preliminary  deci- 
sions. Thus,  the  eventual  decision  may  require  a significant  amount  of  data 
processing. 

False  Alarm  Control : It  is  always  possible  to  improve  the  probability  of 
correctly  deciding  that  an  intruder  is  present  at  the  expense  of  increasing 
the  probability  of  false  alarm.  Since  the  average  level  of  the  signals  out 
of  the  array  may  change  over  a wide  range,  the  appropriate  level  at  which  a 
decision  of  intruder  present  should  be  made  will  also  change.  Thus,  in  ord- 
er to  maintain  the  probability  of  false  alarm  at  a satisfactorily  small 
value,  it  is  nessary  to  adjust  the  decision  levels  as  external  conditions 
change.  This  implies  that  some  computation  must  be  performed  in  order  to 
determine  what  the  appropriate  levels  are  at  any  instant  of  time. 
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Automatic  Seif-Checking:  There  is  a finite  probability  that  any  given  cell 
in  the  array  may  fail.  Such  a failure  may  result  in  a false  alarm,  if  it 
occurs  rapidly,  but  a slow  deterioration  might  never  be  observed  in  this 
way.  Therefore,  it  is  essential  that  some  provision  be  made  for  automati- 
cally checking  the  quality  of  each  cell  on  a periodic  basis.  This  check  re- 
quires some  additional  processing  capability  in  the  system. 

2. 2.£  Data  Storage  Requirements 

The  data  storage  requirement  of  the  processing  system  is  addressed 
separately  because  it  is  a vital  factor  in  determining  what  is  feasible  and 
what  is  not.  It  is  not  intended  to  discuss  the  requirements  in  detail  here, 
but  only  to  outline  the  scope  of  the  problem. 

First  it  may  be  noted  that  to  store  the  data  from  one  stare  time  re- 
quires 8 x 256  = 2048  bits  of  storage.  Next,  suppose  the  reference  signal  ^ 

were  constructed  from  an  average  over  all  of  the  data  for  the  previous  60 
seconds,  and  that  the  stare  time  is  0.1  seconds.  The  total  storage  require- 
ment for  this  one  function  for  this  single  array  would  be 

B = ^ x 2048  = 1.2288  x 106  bits 

It  is  clear  that  some  compromise  must  be  made  between  requirements  that  are 
desirable  and  those  that  are  feasible. 

Other  parts  of  the  system  also  require  data  storage.  For  example,  most 
targets  will  be  in  the  field  of  view  for  several  stare  times.  If  the  target 
is  one  that  is  difficult  to  detect  it  may  be  desirable  to  utilize  the  data 
from  several  stare  times  in  order  to  make  the  decision.  Thus,  it  may  be 
necessary  to  store  data  before  making  the  comparison  and  again  before  making 
the  decision. 
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2. Other  Considerations 

The  intent  of  this  section  is  to  list  other  aspects  of  system  perfor- 
mance that  need  to  be  considered  in  proposing  a tentative  system  design. 

These  aspects  will  not  be  discussed  in  detail;  in  fact,  the  requirements  are 
self-evident  in  most  cases. 

a)  The  system  must  operate  over  a wide  range  of  temperatures.  In 
this  regard,  the  use  of  digital  processing  is  a considerable 
advantage. 

b)  No  manual  adjustments  should  be  required  after  the  initial  in-  ' 

stal lation. 

c)  The  system  must  respond  to  both  slowly  moving  and  rapidly  mov- 
ing intruders.  The  slowest  speed  of  interest  is  on  the  order 

> 

of  .025  m/s  and  greatest  speed  of  interest  is  on  the  order  of 
40  m/s. 

d)  The  system  should  function  under  adverse  weather  conditions 
such  as  fog,  rain,  and  snow. 

I 

e)  The  system  should  be  capable  of  distinguishing  between  small 
intruders,  such  as  rabbits  and  squirrels,  and  larger  intruders 
such  as  dogs  and  men. 

f)  The  system  should  ignore  natural  phenomena  such  as  clouds,  sun 
glint,  lighting,  etc. 

g)  The  system  should  ignore  lights  in  the  field  of  view  that  are 
turned  on  or  off  and  reflections  from  headlights  of  vehicles 
passing  outside  the  field  of  view. 
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h)  System  cost,  reliability,  and  maintainability  are  also  impor- 
tant characteristics,  but  they  cannot  be  adequately  addressed 
in  this  preliminary  study. 

2.2.4  Survey  of  Candidate  Systems 

One  way  of  classifying  different  processing  systems  is  to  base  the 

4 

classification  on  the  form  of  the  reference  and  type  of  comparison  that  is 
made.  On  this  basis  the  most  promising  candidates  for  consideration  may  be 
tabulated  as  shown  in  Table  I.  Any  combination  of  one  reference  signal  and 
one  comparison  method  constitutes  a candidate  system;  thus  defining  eight 
systems . 

Table  _I . Classification  of  Candidate  Systems 
Reference  Signal  Comparison  Method 

1.  Previous  output,  same  cell  A.  Difference 

2.  Adjacent  cell,  same  time  B.  Correlation 

3.  Previous  output  and  adjacent  cell 

4.  Average  of  previous  outputs, 
same  cel l 

It  is  not  the  purpose  of  this  section  to  present  analyses  of  all  candi- 
date systems.  Instead,  some  problems  with  specific  techniques  that  are 
deemed  to  be  sufficiently  serious  to  discourage  further  consideration  will 
be  pointed  out. 

The  difficulty  with  using  the  previous  output  from  each  cell  as  the 
reference  (1)  is  that  for  a slowly  moving  target  the  change  in  cell  output 
from  one  stare  time  to  the  next  may  be  insufficient  to  provide  detection. 
This  difficulty  can  be  alleviated  to  some  extent  by  using  the  output  that 
occurred  n stare  times  earlier  as  the  reference.  A further  difficulty  with 
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either  approach,  however,  is  that  the  reference  is  just  as  noisy  as  the  raw 
data  and,  hence,  larger  signals  are  required  to  achieve  a good  probability 
of  detection  along  with  a small  probability  of  false  alarm. 

When  the  reference  is  derived  from  the  adjacent  cell  (or  cells)  (2) 
there  are  two  difficulties.  One  problem  arises  if  one  cell  is  constantly 
viewing  a discrete  object  that  is  different  from  the  rest  of  the  background 
(e.g.,  a lamppost  or  signpost).  If  no  past  history  is  used,  this  cell  will 
always  declare  a target  present.  The  second  difficulty  is  the  noise  in  the 
reference  as  discussed  above. 

Using  (1)  and  (2)  in  combination  (3)  represents  a more  viable  approach 
than  either  one  above.  However,  it  appears  to  offer  no  advantage  over  (4) 
and,  when  employed  to  extent  that  appears  to  be  necessary,  is  even  more  com- 
plex to  implement. 

The  correlation  method  of  comparison  is  considered  next.  A possible 
method  of  accomplishing  this  is  shown  in  Fig.  2-1.  Although  the  method 
shown  here  involes  correlation  with  the  preceding  sample  only,  it  can  be  ex- 
tended to  include  more  than  one  previous  sample.  The  input,  xk(t),  is  the 
output  of  the 


Fig.  2-1.  Correlation  method  of  comparison. 

kth  sensor  cell  at  the  tth  stare  time,  and  the  reference,  mk,  is  the  mean 
value  of  this  output.  If  a change  of  Ax  occurs  in  the  signal  at  the  f.th 
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stare  time  and  if  there  are  independent  noises  at  the  two  stare  times,  then 


the  output  yk(t)  is 


y^CO  = Cm^+Ax+n^ ( J.)] Cmk+nk ( t-1 )]  - m^ 


= AxCmt+nta-1)]  + mkCnkU)+nkU-1)]  + nL  U)nL  U-1 ) (2-1) 


k k 


The  mean  value  of  this  output  is 


E[yk(H)]  = mkAx 


and  its  variance  is 


(2-2) 


VarCy.  (8.)]  = + (Ax+2m)^  o^ 

k n n 


(2-3) 


where  on  is  the  variance  of  the  noise.  Thus,  the  signal -to-noise  ratio  at 
the  output  is 


(SNR)  = 


(mk Ax) 


a +(Ax+2m.)  oo  _ _ 

n k n n+(Ax+2)2o2 

2 n 
mk  k 


(2-4) 


The  difference  method  of  comparison  is  shown  in  Fig.  2-2. 


x.U)  +• 


ykU) 


Fig.  2-2.  Difference  method  of  comparison. 


Under  the  same  circumstances  as  above,  the  output  in  this  case  is 


(2-5) 


y k C A)  = Cmk  + Ax  + nk(Jl)D  - mk  = Ax  + n^U) 
The  mean  value  of  this  output  is 


E Cy  k C Jt)  3 = Ax  (2-6) 

and  its  variance  is 

VarCy. (1)3  = a2  (2-7) 

k n 


yielding  a signal -to-noise  ratio  of 


(SNR) 

d 


(Ax)2 

2 

a 

n 


(2-8) 


Comparison  of  (2-8)  with  (2-4)  reveals  that  the  difference  method  of 
comparison  yields  a higher  signal -to-noise  ratio  under  all  circumstances. 
This  aspect,  combined  with  the  fact  that  the  correlation  method  is  more  dif- 
ficult to  implement,  appears  to  dictate  the  use  of  the  difference  method. 
Thus,  the  candidate  system  that  appears  to  have  the  greatest  promise  is  4A. 
It  is  this  system  that  is  analyzed  in  more  detail.. 
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2.3  PROPOSED  SYSTEM  DESIGN 


2. 3^  1^  Genera l Descript  ion 

The  system  that  is  proposed  for  analysis  and  evaluation  is  one  in  which 
the  reference  signal  is  constructed  by  averaging  the  output  of  each  cell 
over  a time  interval  that  is  long  compared  to  the  time  that  an  intruder  is 
likely  to  remain  in  the  field  of  view,  and  the  comparison  is  made  by  sub- 
tracting the  reference  signal  from  the  most  recent  output  of  each  cell.  A 
block  diagram  of  this  system  is  shown  in  Fig.  3-1,  and  the  operation  per- 
formed by  each  block  are  discussed  in  the  following  paragraphs. 

Sensor  Array.  This  block  contains  the  IR  detector  cells,  the  CCD,  and  the 
associated  circuitry  required  to  read  out  the  signals  from  the  CCD.  The 
stare  time  adjustment  is  also  accomplished  here. 

A/D  Converter.  Converts  each  output  from  the  CCD  into  an  8-bit  binary  se- 
quence. Since  there  are  4u  seconds  for  each  converson,  the  bit  rate  out  of 
the  A/D  converter  is  2Mbps  during  read  out.  Data  is  available  for  1.024  ms 
during  each  stare  time.  The  output  of  the  A/D  converter  is  denoted  as  x^i), 
which  is  a vector  having  256  elements  (the  number  of  cells  in  the  array), 
each  element  of  which  is  an  8-bit  word. 

Background  Spatial  Average.  The  outputs  from  a subset  of  cells  that  are 
more  or  less  uniformly  spaced  in  the  array  are  averaged  as  a means  of  ob- 
taining a signal  that  can  be  used  to  adjust  the  stare  time.  If  D such  cells 
are  used,  they  can  be  defined  by  a vector 

uj  = (0,1, 0“*  0,1,0) 

in  which  D of  the  elements  are  unity  and  the  rest  are  zero.  The  average  at 
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Sensor  I Stare  Time 


gure  3*1.  Block  diagram  of  proposed  system. 


each  stare  time  is  given  by 


c(l)  = i uj  (3-1) 

Stare  Time  Control  . The  value  of  cU)  is  used  to  make  adjustments  in  the 
stare  time.  A possible  rule  for  doing  this  is  as  follows: 

1)  Set  two  thresholds,  and  where  c2  is  somewhat  lower  than  the 
level  at  which  saturation  occurs  and  c^  is  on  the  order  of  250  mv 
less  than  c^. 

2)  Adjustments  to  stare  time  are: 

C1  £ c(D  < c^:  No  change 
c(t)  > c2:  Reduce  stare  time  by  2 
cU)  < c^:  Increase  stare  time  by  2 
The  value  of  D must  be  large  enough  that  the  presence  of  a target  in 

one  of  the  cells  entering  into  the  average  of  (3-1)  will  not  significantly 

affect  the  average.  It  is  believed  that  D = 16  is  sufficient.  In  this  case 
the  elements  of  u^  would  be 

uk  3 1 , k 3 8 ♦ 16j  , j 3 0,  1,  2,  ...,  15 

3 0 , otherwise  . 

It  may  be  noted  that  a decision  to  change  or  not  change  the  stare  time 
is  made  after  each  stare  time  and  that  there  is  no  long-  term  time  average. 
Such  a time  average  is  not  necessary  because  the  hysteresis  in  the  decision 
rule  prevents  oscillation  across  the  boundary. 

Target  Time  Average.  The  possibility  of  averaging  the  target  signal  over 
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several  stare  times  is  included  in  order  to  take  advantage  of  the  fact  that 
most  targets  remain  in  the  field  of  view  for  several  stare  times.  Subse- 
quent evaluation  of  the  system  performance  indicates  that  such  averaging  is 
probably  not  necessary  under  most  circumstances,  but  the  possibility  is  re- 
tained here  in  order  to  insure  the  greatest  system  flexibility. 

The  target  time  average  for  each  cell  is  the  average  of  the  M most  re- 
cent outputs  for  that  cell.  Thus,  it  may  be  defined  as 

. M-1 

yU)  = ^ (3-2) 

m=0 

If  this  operation  is  not  used,  then  11=1  and  y_(l)  = x_(l).  If  this  operation 
is  used,  the  value  of  M should  be  no  larger  than  is  necessary  to  span  the 
number  of  stare  time  that  the  fastest  moving  small  target  is  in  the  field  of 
view. 

Background  Time  Average.  It  is  necessary  to  generate  a reference  signal  for 
each  cell.  The  method  employed  here  is  to  perform  a long-time  average  on 
the  background  signal  generated  by  each  cell.  Since  the  important  criterion 
is  the  length  of  time  over  which  the  average  is  taken  rather  than  the  number 
of  samples  actually  averaged,  the  proposed  method  performs  the  average  over 
a set  of  stare  times  separated  in  time  by  P stare  time  intervals  in  order  to 
reduce  data  storage  requirements.  Thus,  the  background  time  average  can  be 
expressed  as 

1 N 

m(i')  = jj-  £ jjCPU'-nn  (3-3) 

n=1 

in  which  the  index  l*  * l/P  corresponds  to  every  Pth  stare  time.  The  total 
time  interval  over  which  this  average  is  taken  is  PNt$  and  should  be  long 
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compared  to  the  time  that  a target  remains  in  the  field  of  view.  This  is 
because  any  target  that  remains  in  the  field  of  view  this  long  becomes  part 
of  the  background  average  and  is  no  longer  detectable. 

Although  the  product  PN  is  the  important  parameter  for  this  operation, 
the  parameter  N does  play  a role  in  determining  the  variance  of  the  average. 
In  order  to  reduce  this  variance  to  a value  significantly  smaller  than  the 
variance  of  the  target  signal,  it  is  necessary  that  N»M. 

Store  and  Re-sample.  Since  a new  reference  signal  is  created  only  once 
every  P stare  times,  but  a reference  is  needed  every  stare  time,  it  is 
necessary  to  store  the  reference  signal  and  re-sample  it  at  the  higher  rate. 
Thus,  there  will  be  P successive  vectors,  mU),  that  are  identical. 

Comparison  Operation.  The  comparison  operation  in  this  sytem  is  simple  vec- 
tor subtraction.  Thus, 

£(l)  = y(l)  - m(l)  (3-4) 

in  which  each  element  of  z_(.i)  is  the  difference  between  the  output  of  that 
particular  cell  and  the  background  average  for  that  particular  cell.  In  the 
absence  of  a target,  each  element  of  £(l)  is  an  8-bit  word  representing  a 
zero-mean  random  variable. 

Absolute  Value  Operation.  The  change  in  the  mean  value  of  £(l)  that  results 
from  a target  entering  the  field  of  view  can  be  either  positive  or  negative. 
In  order  to  avoid  the  need  Tor  establishing  two  thresholds,  it  convenient  to 
perform  an  absolute-value  operation.  Thus,  wU)  is  a vector,  each  element 
of  which  is  the  absolute  value  of  the  corresponding  element  of  £(i).  That 
is,  the  kth  element  is 
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. 


i. 


wRU)  = |zkU)| 


(3-5) 


Thresho 1 d Computation . The  decision  as  to  the  presence  or  absence  of  a tar- 
get is  based  on  comparing  each  element  of  w(l)  with  a threshold  bU).  The 
proper  value  for  the  threshold  depends  upon  the  variance  of  w^l).  As  is 
shown  subsequently,  this  variance  is  uniquely  related  to  the  mean  value  of 
£(*.).  Thus,  this  vector  provides  the  input  necessary  to  calculate  the 
correct  threshold  values.  The  threshold  also  depends  upon  the  desired  pro- 
bability of  false  alarm.  The  method  of  calculation  is  described  in  a subse- 
quent section. 

Threshold  Operation.  This  operation  compares  each  element  of  w^(t)  with  the 
corresponding  element  of  the  threshold  vector,  t^OO.  If  the  threshold  is 
exceeded,  a logical  1 is  produced,  if  not  a logical  0 is  produced.  These 
results  are  used  to  form  a vector  \/(t),  each  element  of  which  is  a 3-bit 
word  corresponding  to  the  outcomes  of  the  threshold  comparisons  for  three 
adjacent  cells.  If  the  result  of  the  ith  comparison  is  a^U)=  0 or  1,  then 
the  vector  v(t)  is 


a^i) 

32(0 

a3(t)- 

a^t) 

a^U) 

a.  U) 

H 

3j(l) 

• 

a4U) 

• 

a^U) 

• 

• 

a253U) 

• 

a254U) 

» 

a255U) 

a 254  ( 4 ) 

3 255 

a256U) 

This  will  be  a (254x1)  vector  unless  the  extremal  cells  of  adjacent  arrays 
are  incorporated  in  this  operation,  in  which  case  it  would  be  (256x1).  This 
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latter  mode  of  operation  implies  a central  processor  operating  on  the  out- 
puts of  all  arrays  rather  than  a microprocessor  associated  with  each  array. 


Decision.  In  the  absence  of  a target,  most  of  the  elements  of  v^l)  are 
zeros.  When  a detectable  target  enters  the  field  of  view,  the  elements  of 
vU)  corresponding  to  the  appropriate  cells  will  change  to  ones  and  will 
remain  ones  for  as  many  stare  times  as  the  target  is  in  the  field  of  view. 
In  order  to  provide  discrimination  between  desired  and  undesired  targets, 
the  eventual  decision  is  based  on  the  states  of  n/3  successive  threshold 
comparisons  in  each  of  the  three  cells  associated  with  any  element  of  v/O. 
A decision  of  "target  present"  is  made  of  t or  more  of  the  comparison  in  any 
block  of  n are  one. 

As  an  illustration  suppose  that  n=12  and  t=8.  The  block  of  digits 
corresponding  to  cells  i,  (i+1)  and  (1+2)  after  the  Hh  stare  time  would  be 

a.U) 
a^ (l-l) 
a-U-2) 
a.U-Z) 

If  8 or  more  of  these  digits  are  ones,  a decision  of  "target  present"  is 
made.  If  not,  a decision  of  "no  target"  is  made.  This  decision  operation 
is  performed  once  every  stare  time. 

Self-Checking  Feature.  An  essential  part  of  the  proposed  system  is  the 
ability  to  periodically  check  its  own  operation  and  give  an  indication  that 
all  cells  are  functioning  and  that  all  of  the  computation  is  being  performed 
correctly.  This  feature  is  indicated  on  the  block  diagram  of  Fig.  3-1  by 
the  input  to  the  decision  block  labeled  "Initiate  Zero  Count",  and  function 


v.U)  = 


ai+1(l) 

ai+2U) 

i+1u-1> 

ai+2(i-1) 

i+1 (i-2) 

ai+2U-2) 

i+1 (t'3) 

ai+2U“3) 

is  performed  as  follows: 

1)  At  appropriate  time  intervals  (e.g.,  once  every  five  minutes)  a 
heater  near  the  sensor  array  is  activated.  This  heater  should 
have  a short  time  constant  (<0.1  second)  and  be  positioned  so  that 
all  cells  of  the  sensor  are  illuminated  equally.  The  temperature 
differential  established  by  the  heater  should  be  slightly  greater 
than  that  necessary  to  result  in  a decision  of  target  present. 

2)  At  the  time  the  heater  is  activated,  the  decision  logic  should  be 

altered  to  count  the  number  of  zeros  in  every  set  of  n observation 

rather  than  the  number  of  ones,  and  t is  set  equal  to  n/3. 

3)  If  the  system  is  operating  properly,  every  threshold  should  be  ex- 
ceeded and  there  should  be  no  zeros.  If  any  cell  is  not  function- 

ing, it  will  produce  zeros  and  an  alarm  will  be  sounded  for  that 
cell.  Likewise  if  the  thresholds  or  the  computations  are  in- 
correct, and  if  this  error  has  not  previously  caused  an  alarm,  it 
will  do  so  when  the  decision  is  complemented. 

Anal  ysis  of  System  Operation 

Although  the  system  operation  described  above  appears  to  be  complex, 
because  of  the  large  amount  of  data  that  needs  to  be  handled  in  each  step, 
the  analysis  of  this  operation  is  really  quite  straightforward.  Because  the 
operations  performed  on  the  output  of  each  cell  of  the  sensor  are  the  same, 
the  analysis  is  carried  out  for  a single  cell.  To  further  simplify  the 
analysis,  the  output  of  each  cell  is  considered  to  be  a continuous  random 
variable  rather  than  a quantized  one.  Hence,  the  analysis  is  for  an  analog 
system  rather  than  a digital  one.  Because  the  quantizing  error  is  much 
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smaller  than  the  expected  signal  levels,  either  with  or  without  a target, 
this  approximation  should  be  a good  one.  Furthermore,  it  should  be  em- 
phasized that  the  computer  evaluation  of  the  system  did  perform  the  ap- 
propriate quantization  so  that  even  this  approximation  did  not  exist  in  that 
study. 

When  there  is  no  target  present,  the  output  of  any  one  cell  in  any  one 
stare  time  is  a random  variable,  Xq-  Because  of  the  Poisson  distribution  of 
the  number  of  electrons  contributing  to  this  output,  the  random  variable  Xq 
has  a mean  value  of 


= eNB  (3-8) 

where  N is  the  average  number  of  electrons  as  defined  by  (1-11)  or  (1-12), 

D 

and 


(3-9) 


-1 9 

where  e = 1.6x10  coulombs  is  the  electronic  charge  and  Cg  is  the 
equivalent  capacitance  of  the  CCD  transfer  circuit  that  determines  the  out- 
put voltage.  For  example,  if  Cg  = 0.5  pF,  as  assumed  in  our  calculations, 
then  B = 3.2x10-7.  Because  the  variance  of  a Poisson  distribution  is  the 
same  as  its  mean,  the  random  variable  Xq  has  a variance  of 


Because  N0  is  so  large,  it  is  reasonable  to  assume  that  xq  is  a con- 
tinuously distributed  random  variable  with  a Gaussian  probability  density 
function.  It  is  further  assumed  that  the  variations  in  Xq  are  independent 
from  one  stare  time  to  the  next  because  the  variations  in  the  number  of 
electrons  are  independent.  Thus,  it  is  possible  to  express  the  means  and 


i 
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variances  of  the  target  time  average  and  the  backgrond  time  average,  when  no 

2 

target  is  present,  in  terms  of  mn  and  on.  The  target  time  average  is 


, H-1 

'■i  5 


(3-10) 


m=0 


where  x is  the  cell  output  at  the  mth  stare  time  preceding  the  current  one. 
m 

It  follows  that  the  mean  value  of  y is 


. M-1  r . M-1 

y0  = ECy;i  =F  £.  E[xm]  = tt  ^ 
m=0  L J m=0 


m0  = m0 


(3-11) 


Similarly,  because  of  independence,  the  variance  of  y is 


M-1  M-1 


yo 


E[(y“/)?]=^  ?n  EC<Vx0)(V-x0)] 

*-  J M m=0  m ' =0 


1 


M-1 


77  £ EC<Vx0)2:1  = "7  £ 

M m=0 


M-!  2 oQ 

7 m=0  °°  = 


(3-12) 


The  background  time  average  is 


1 ™ 

m = i L xn 

n=1 


(3-13) 


where  xn  is  the  cell  output  at  nPth  stare  time  preceding  the  current  one. 
By  analogy  to  (3-11)  and  (3-12),  the  mean  and  variance  of  this  average  is 


m0  = ""O 


(3-14) 
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(3-15) 


The  output  of  the  comparison  operation  is 

z = y-m  (3-16) 

If  it  is  assumed  that  P>M,  a condition  that  will  almost  certainly  be  true, 
then  Y and  M contain  no  outputs  from  the  same  stare  times  and,  hence,  will 
be  independent.  Thus,  the  mean  and  variance  of  z,  when  there  is  no  target, 
are 


= y -m„  = m -in  = 0 
o oo  oo 


(3-17) 


and 


2 2 2 2 

a = o +a  = a 
zo  yo  mo  o 


[ff  + tt] 


(3-18) 


Since  z is  also  a Gaussian  random  variable,  the  probability  that  its 
magnitude 


w = |z| 


(3-19) 


exceeds  a threshold  value,  b,  is  simply 


pQ  = PrCz>b]  + Pr[z<-b]  = 2 Q 


zo 


(3-20) 


where 


Q(c)  = — — f e"u2/2du 

V7n  a 

is  the  complementary  Gaussian  distribution  function.  The  relation 
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Q (a)  = l-Q(-a) 


(3-21) 


is  also  used  in  evaluating  that  portion  of  (3-20)  that  pertains  to  z<-b.  the 
probability  pq  is  related  to  the  probability  of  false  alarm  and  this  rela- 
tionship is  derived  subsequently. 

When  a target  enters  the  field  of  view  of  any  cell,  both  the  mean  and 
variance  of  the  cell  output  change.  Specifically,  the  new  mean  value  is 


m1  = 8Nt  (3-22) 

where  is  the  average  number  of  electrons  as  defined  by  (1-15)  or  (1-16). 
Similarly,  the  new  variance  becomes 


o 


2 

1 


(3-23) 


However,  these  values  will  change  from  stare  time  to  stare  time  as  the  tar- 
get moves  into  and  out  of  the  field  of  view;  that  is,  the  fractional  area  6 
is  a function  of  time.  It  is  necessary,  therefore,  to  expand  t • above  no- 
tation to  let  mf  be  the  mean  value  of  the  cell  output  at  the  nth  stare  time, 
• 2 

when  a target  is  present,  and  is  the  variance  at  the  same  time. 

The  target  time  average  at  the  nth  stare  time,  by  analogy  to  (3-11),  is 


M-1 

E 


m=0 


n-m 


(3-24) 


and  the  variance,  by  analogy  to  (3-12),  is 
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n-m 


(3-25) 


The  background  time  average  can  be  handled  in  a simitar  fashion  to 


yield  a mean  and  variance  of 


(3-26) 
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(3-27) 


Again  under  the  assumption  that  P > M,  the  mean  and  variance  of  z , the 
output  of  the  comparison  circuit,  are 


1 M~.1  1 N, 
m - -jj  / m • tt  V m _ 
ZH  M i-m  N fc-nP 

m=U  n=l 


and 


(3-28) 
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(3-29) 


It  may  be  noted  that  the  mean  value  is  no  longer  zero.  Hence,  the  probabil- 
ity that  w = I z | will  exceed  the  threshold  b is  not  the  same  as  it  was 
when  there  was  no  target  present.  The  situation  is  represented  graphically 


Fig.  3-2.  Illustrating  the  probability  of  exceeding  the  threshold. 
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the  threshold  is 


pi  ■ pr[vb]  * Pr[v-b] 


= Q 


fb-m 

m +b 
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+ Q 

N 

1 
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= 1-Q 
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+ Q 
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(3-30) 


The  probability  is  related  to  the  probability  of  detection  and  this  rela- 
tionship is  presented  subsequently. 

It  may  also  be  noted  that  (3-30)  includes  the  case  of  no  target  simply 

by  setting  m =0  and  Y = v 

Zl  xz t xzo 

For  computational  purposes  it  is  convenient  to  calculate  the  signat- 
to-noise  ratio  for  each  stare  time  since  this  parameter  is  the  one  that  is 
most  significant  in  determining  the  performance.  The  signal -to-noise  ratio 
for  the  j,th  stare  time  is 


(m  - m ) 

(SNR)  = — & — r-2- 
l 2 


(3-31) 


From  this  definition,  an  average  signal -to-noise  ratio  for  the  target  time 
average  can  be  expressed  symbolically  as 


i i M_1  r 

V<s"Rw  V 


(SNR) 


m=0 


j-m 


(3-32) 


Similarly,  an  average  signal-to-noise  ratio  for  the  background  time  average 
can  be  expressed  symbolically  as 

34 


-- 


= ff  £ V<SN"VnP 

Also  let 


(3-33) 


(3-34) 


and  define 


M-1 

ST  = 'T  ^ 
m=0 
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(3-35) 
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Final ly,  define 


and 


d = Q_1(po/2) 


(3-36) 


(3-37) 


where  pQ  is  the  desired  value  of  pQ  as  determined  by  the  desired  probability 
of  false  alarm  and  Q (•)  is  the  inverse  Q-function.  Thus,  from  (3-20) 
threshold  value,  b,  is 


do 


•V 


1 + 1 
W TT 


(3-38) 
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If  all  of  the  above  definitions  are  substituted  into  (3-30),  the  result 


is 


pt  = 1-Q 


VC 

s 

(SNR)t  - ^ 

(SNR)b 

-d 

1 

G(St+sb)+1 

+ Q 


VC 


\ 


^T(SNR)t  - ^f(SNR)E 
■^G(ST+SB)+1 


+d 


(3-39) 


The  advantage  to  this  formulation  can  be  seen  by  writing  the  signal -to-noise 
ratio  at  the  Ath  stare  time,  (3-31),  as 


(SNR)  = 
A 


mrmo 


bNta~bNb 

b2n„ 


nta~nb 


(3-40) 


and  the  parameter  S^,  (3-34),  as 


r _ mA_mo  BNTi~BNB  NTa“NB 
bA  m ~ BN 


B 


N 


B 


(3-41) 


It  is  noted  that  these  parameters  can  be  expressed  entirely  in  terms  of  the 
number  of  electrons  and  do  not  require  a knowledge  of  the  parameter  B which 
converts  electrons  to  volts.  Thus,  (3-39)  is  general  and  would  be  applica- 
ble to  any  array  of  sensors  having  any  transfer  efficieny  for  the  CCD. 

It  is  true  that  a knowledge  of  B is  necessary  to  actually  evaluate  the 
threshold  b as  defined  by  (3-38).  However,  the  value  of  b is  only  required 
in  the  operating  system  and  is  not  required  in  calculating  p^.  In  the 

operating  system,  o (which  is  a factor  of  b)  is  proportional  to  J m and  an 
o 1 o 

estimate  of  mQ  is  continuously  available  at  the  output  of  the  background 
time  average.  Thus,  b is  proportional  to  ^m^  and  the  constant  of  propor- 
tionality is  obtained  from  actual  adjustment  in  the  operating  system  to 
achieve  the  desired  performance. 


36 


r 


- . . _ 


and 


The  final  step  in  the  analysis  of  system  operation  is  to  relate  p 

I 

p^  to  the  probability  of  false  alarm  and  the  probability  of  detection.  As 

discussed  in  Sec.  2.3.1/  the  final  decision  is  based  on  observing  the  number 

of  times  the  threshold  has  been  exceeded  in  the  n/3  most  recent  stare  times 

in  3 adjacent  channels.  If  t or  more  of  these  observations  exceed  the 

threshold,  a decision  of  "target  present"  is  made. 

For  the  case  in  which  a target  is  not  actually  present/  the  probability 

of  false  alarm  is  related  to  p by  the  binomial  distribution  since  p is  the 

o o 

probability  that  noise  alone  will  exceed  the  threshold.  Thus,  the  probabil- 
ity of  false  alarm  is  given  by 


P 


F 


(3-42) 


The  importance  of  this  relation  lies  in  the  fact  that  extremely  small  values 
of  Pp  can  be  obtained  with  only  moderately  small  values  of  pq. 

By  way  of  illustrating  the  above  point,  suppose  that  a 5 km  perimeter 
is  viewed  by  130  arrays  each  containing  256  sensors  and  each  sensor  views 
0.15  m of  the  perimeter.  Thus,  there  are  a total  of  33280  sensors.  If  only 
one  false  alarm  per  day,  due  to  system  noise,  is  desired,  and  if  the  stare 
time  is  0.1  second,  then  the  probability  of  false  alarm  for  each  sensor 
would  have  to  be 

d _ 0.1  -x/.o^in-ll 

PF  - 33,280(741(76(10)  3-48  x 1U 

which  is  an  extremely  small  value.  However,  by  selecting  n=12  and  t=8,  the 
PQ  required  to  achieve  this  is  0.0229. 

When  a target  enters  the  field  of  view,  not  all  of  the  p^  will  have  the 
same  value  since  they  correspond  to  different  observation  times  and  d i f — 
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ferent  cells.  Thus,  there  will  be  a set  of  n values  of  p . Designate  these 
n values  as  and  let  = 1-p^.  The  probability  of  detection  is  then 

given  by 

t-1  j 

PD  = 1 “ £ (pi  *•••  P|<  Pf  PQ>  (3-43) 

j=0  9 

where  the  second  sum  is  over  all  combinations  of  p..,  q..  containing  j value 
of  p.j  and  n-j  values  q^.  To  illustrate  this  more  clearly,  let  n=3.  Then 

j=o 

£ (Pi  '"/  Pf  ••*>  = qi  P2  P3 

j=1 

Z ( } = P1  q2  q3  + p2  q1  q3  + p3  P1  ^2 

J 

j=2 

Z < ) = P1  P2  q3  + P1  P3  q2  + P2  P3  P1 

jf3  , , 

Z < > - p,  p2  p3 

In  addition,  if  t=2  in  this  example,  then  j=0  and  j=1  appear  in  (3-43)  and 
the  probability  of  detection  becomes 

PD  * 1 • fl1  "2  q3  - P1  q2  q3  - p2  q1  q3  - p3  q1  q2 

Clearly,  this  becomes  a very  involved  computation  for  large  values  of  n and 
t. 

Again,  the  importance  of  this  result  is  that  values  of  very  close  to 
1 can  be  obtained  with  much  smaller  values  of  p . To  illustrate  this  point 
further,  suppose  it  is  desired  to  achieve  a probability  of  detection  of 
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0.999  and  all  of  the  values  are  the  same  (e.g.  a target  has  completely 
covered  the  field  of  view  of  3 cells  for  at  least  n/3  stare  times).  Then 
for  n=12  and  t=8  again,  the  required  value  of  p^  is  only  0.784. 

_2.3_.  3.  Comparison  with  An  Optimum  System 

It  is  of  interest  to  compare  the  performance  of  the  proposed  system 
with  that  of  an  optimum  system.  It  is  to  be  expected  that  the  optimum  sys- 
tem will  perform  somewhat  better  for  at  least  three  reasons: 

1)  The  optimum  detector  requires  nonlinear  processing  of  the  observed 
data,  while  the  proposed  system  is  linear  except  for  the  absolute 
value  operation  on  z. 

2)  The  optimum  detector  requires  that  data  from  all  cells  and  for  all 
observation  times  be  processed  to  make  a single  decision.  The 
proposed  system  uses  data  from  only  3 cells  and  M observation 
times  for  each  decision. 

3)  The  thresholds  of  the  optimum  are  assumed  to  be  known  exactly  from 
a complete  knowledge  of  the  background  and  target  characteristics. 
In  the  proposed  system  the  threshold  is  computed  from  an  estimate 
of  the  variance  of  the  observed  processes. 

The  analysis  of  the  optimum  detector  is  presented  in  Appendix  A in  its 
general  form.  Because  of  computational  difficulties  in  integrating  multidi- 
mensional Gaussian  density  functions,  the  only  situation  for  which  it  is 
feasible  to  make  a numerical  comparison  is  the  single  observation  case.  The 
special  result  for  the  optimum  system  for  this  condition  is  also  given  in 
Appendix  A.  For  the  proposed  system,  this  condition  corresponds  to  setting 
M=n=t=1.  The  resulting  comparison  is  shown  in  Fig.  3-3  by  displaying  the 
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probability  of  detection  as  a function  of  (SNR)^  for  a false  alarm  proba- 
bility of  6.94  x 10  It  is  evident  that  the  proposed  system  is  subop- 

timum, but  only  by  a small  amount  x - less  than  1 dB. 

It  is  not  known  whether  the  proposed  system  is  better  or  worse  relative 
to  the  optimum  system  when  more  than  one  observation  is  considered,  because 
the  optimum  system  cannot  be  computed  for  this  case  without  utilizing  an  ex- 
cessive amount  of  computing  time.  Furthermore,  since  the  numerical  evalua- 
tion of  the  proposed  system  performance  produced  results  that  exceeded  ex- 
pectations, the  expenditure  of  such  computing  effort  for  the  sake  of  making 
further  comparisons  does  not  seem  justified. 

2^. 3.4.  Selection  of  System  Parameters 

In  order  to  evaluate  the  performance  of  the  proposed  system  it  is  first 
necessary  to  select  the  operating  parameters.  This  was  accomplished  by 
varying  the  parameters,  one  at  a time,  and  computing  the  resulting  probabil- 
ity of  false  alarm  and  probability  of  detection  for  two  "worst  case"  tar- 
gets. These  worst  case  targets  are  designated  as 

I.  The  most  detectable  target  for  which  a decision  of  "no  target"  is 
desired.  This  turned  out  to  be  the  target  model  for  a squirrel 
against  the  most  contrasting  background. 

II.  The  l east  detectable  target  for  which  a decision  of  "target"  is 
desired.  This  is  the  target  model  for  a crawling  man  against  the 
least  contrasting  background. 

The  parameter  variations  for  each  of  these  targets  are  described  in  the  fol- 
lowing paragraphs. 
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Figure  3-3.  Single  Observation  Probability  of  Detection 


Number  of  Target  Samples/  M.  The  parameter  M was  varied  from  1 to  10.  It 
was  found  that  M=1  would  produce  satisfactory  results  for  both  worst  case 
targets,  and  that  larger  values  of  M were  undesirable  when  the  target  was 
moving  fast  enough  to  pass  through  the  field  of  view  in  one  or  two  stare 
times.  Hence,  M=1  is  selected  for  use. 

Number  of  Background  Sampl es,  Nu  The  value  of  N was  varied  from  1 to  80. 
It  was  found  that  for  N < 15,  the  background  average  is  seriously  affected 
by  the  presence  of  the  target.  With  values  of  N 40,  the  target  had  little 
effect  on  the  background  average.  In  order  to  minimize  data  storage,  a 
value  of  N=40  is  selected. 

Separation  of  Background  Samples,  P.  the  value  of  P was  tried  at  both  5 and 
10.  It  was  found  that  the  important  parameter  was  really  the  product  NP.. 
Thus,  P=10  and  N=40  produced  the  same  result  as  P=5  and  N=80.  In  order  to 
minimize  data  storage,  P=10  is  selected.  It  should  be  noted  here  that  this 
phase  of  the  study  was  done  before  the  microprocessor  design  was  begun. 
After  consideration  of  the  microprocessor  it  became  apparent  that  N should 
be  a power  of  2.  Thus,  in  the  actual  implementation  of  the  proposed  system 
it  is  recommended  that  N=32  and  P=12  be  used.  Nevertheless,  all  of  the  com- 
putations reported  here  as  based  on  N=40  and  P=10. 

Detection  Threshold,  d^.  Since  both  the  probability  of  false  alarm  and  pro- 
bability of  detection  are  critically  dependent  on  d,  this  parameter  was 
varied  only  in  the  range  from  6 to  7.  It  was  found  that  d=6.168,  which 
corresponds  to  Pp=6.94x10_1®,  provided  satisfactory  values  for  the  probabil- 
ity of  detection. 
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Decision  Threshold,  t^.  The  number  of  observations  used  to  provide  a final 
decision  is  arbitrarily  set  at  n=12,  because  this  is  a multiple  of  3 and  is 
not  too  large  to  prohibit  efficient  computation.  The  threshold,  t,  was  then 
varied  from  6 to  9.  It  was  found  that  in  order  to  detect  target  II,  t must 
be  less  than  9;  and  in  order  to  reject  target  I,  t must  be  greater  than  7. 
Thus,  t=8  is  selected. 

A number  of  parameters  were  not  varied  in  these  tests.  These  are 
described  below. 

Stare  Time,  tg.  A value  of  ts=0.1  second  was  found  to  be  satisfactory  to 
avoid  saturation  and  still  provide  adequate  signal  levels  for  both  target 
cases. 

Depression  Angle,  The  detector  array  was  aimed  at  an  angle  0^  = 5.65° 
below  the  horizon  in  order  to  provide  a range  on  the  order  of  100  meters. 

Quantization  Levels.  Since  8 bits  yields  256  quantization  levels,  and  this 
number  yields  a quantizing  error  on  the  same  order  of  magnitude  as  the  sys- 
tem noise,  this  value  was  selected.  Furthermore,  8 bits  is  a convenient 
number  for  most  microprocessors. 

Optical  Parameter.  The  detector  bandwidth  used  for  these  tests  was  3.4  < X 
< 4.2  um.  The  lense  system  was  assumed  to  have  a focal  length  of  1.7  cm  and 
an  f-number  of  1.  Atmospheric  attenuation  was  assumed  to  be  0.01  dB/m. 
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2.4.  EVALUATION  OF  PROPOSED  SYSTEM 


2.4.1.  Differential  Sens  itivity 

One  way  of  evaluating  the  performance  of  the  proposed  system  is  to 
determine  the  amount  that  a target's  temperature  and/or  emissivity  must 
differ  from  the  background  in  order  for  the  target  to  be  detectable.  For 
small  changes  in  temperature  or  emissivity,  the  change  in  probability  of 
detection  is  expected  to  be  a linear  one.  However,  the  fraction  of  the 
field  of  view  covered  by  target,  6,  also  is  important. 

A convenient  graphical  presentation  of  this  differential  sensitivity  is 
to  plot  lines  of  constant  probability  of  detection  on  a plane  whose  coordi- 
nates are  6Ae  and  6AT,  where 

Ae  = cT  - e8 

and 

AT  = TT  ‘ TB 

The  slopes  and  positions  of  these  lines  of  constant  P^  depend  upon  the  back- 
ground temperature  and  emissivity,  as  well  as  the  parameters  of  the  system. 
Hence,  the  data  from  which  the  plots  are  constructed  must  be  obtained  by 
direct  calculation  of  P^  for  a series  of  small  values  of  Ae  and  AT.  The 
results  of  this  computation  are  displayed  in  Fig.  4-1  through  4-18  for  three 
different  background  temperatures,  three  different  background  emissivities, 
and  for  both  night  and  day. 

Several  points  are  worth  noting  in  connection  with  these  results: 

1)  The  lines  of  constant  P^  have  a negative  slope  for  night  operation 
and  a positive  slope  for  day  operation.  This  is  because  the  solar 
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reflected  signal  is  predominate  in  the  day  and  this  signal  is  pro- 
portional to  1-e.  This  also  suggests  that  at  some  solar  illumina- 
tion level  the  sensitivity  will  be  independent  of  the  temperature 
change . 

2)  In  every  case  there  is  a curve  for  which  PD=0.  This  is  a condi- 
tion in  which  the  change  in  temperature  and  change  in  emissivity 
exactly  compensate  one  another. 

3)  For  night  operation,  the  sensitivity  i ncreases  as  the  background 
emissivity  and  the  background  temperature  increase. 

4)  For  day  operation,  the  sensitivity  increases  as  the  background  em- 
issivity increases  and  the  background  temperature  decreases . 

2. ^2.  Dependence  on  Signal -to-Noise  Ratio 

A second  method  of  evaluating  the  performance  of  the  proposed  system  is 
to  determine  the  overall  probability  of  detection,  P^,  as  a function  of  the 
signal -to-noise  ratio  at  each  observation.  This  relationship  has  been  com- 
puted for  the  parameters  selected  in  Sec.  2.3.4  and  is  displayed  in  Fig. 
4-19.  For  purposes  of  comparison,  the  single-observation  detection  curve  of 
Fig.  3-2  is  also  shown.  It  is  apparent  that  the  use  of  multiple  observa- 
tions significantly  increases  the  steepness  of  the  detection  curve.  This  is 
desirable  because  it  increases  the  probability  of  detection  for  those  tar- 
gets that  should  be  detected  and  reduces  the  probability  of  false  alarm  for 
those  targets  that  should  not  be  detected. 
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Figure  4-2.  Differential  sensitivity  plo 
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Figure  I4-9.  Differential  sensitivity  plot.  Parameter  is  probability  of  detection. 
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Figure  Probability  of  detection  as  a function  of  s ignal -to-no ise  ratio  in  each  observation. 


2.4.3.  Target  and  Background  Model s. 

The  purpose  of  this  section  is  to  tabulate  the  models  for  the  targets 
and  backgrounds  that  were  used  to  evaluate  the  performance  of  the  proposed 
system. 

The  targets  are  modeled  by  rectangular  paral l el epipeds  whose  dimensions 
are  seclected  to  approximate  the  target's  size.  They  are  assumed  to  enter 
the  detector's  field  of  vieu  with  a specified  velocity  and  a specified  angle 
as  shown  in  Fig.  4-20.  Their  emissivity  and  temperature  are  selected  to 
correspond  to  the  physical  characteristics  of  the  surface.  The  data  were 
obtained  from  a variety  of  the  references  listed  at  the  end  of  this  report. 


Fig.  4-20.  Motion  of  target  through  the  field  of  view  of  a single  cell. 
For  some  target  models  it  was  found  to  be  appropriate  to  consider  dif- 
ferent parts  of  the  surface  to  be  at  different  temperatures  or  emissivities. 
These  characteristics,  as  well  as  the  dimensions  and  emissivities  and  speeds 
are  tabulated  in  Table  4-1  for  all  of  the  targets  considered. 

The  field  of  view  that  was  used  for  most  calculations  is  shown  in  Fig. 
4-21  as  a ground  projection.  This  is 
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Table  4-1.  Tabulation  of  Target  Models 


Target  Dimensions  (m) 
Height 
Width 
Length 


Temperature,  Emissivity 
and  Percent  Area 


T,(°K)  e 


Speed 

Range 

m/s 


Squi rrel 


305  (Body) 
288  (Tail) 


Man 

Crawl ing 


Man 

Wal king 


Bl  ac  k 
Car 


305  (Skin)  .98  0.025 

288  (Clothing)!. 80  1.0 


Horse 


Bl  ack 
Truck 


305 

(Sk in) 

288 

(Clothing) 

305 

(Hood) 

288 

(Body) 

305 

(Hood)  . 

288 

(Body) 

305 

305 

(Hood) 

288 

(Body) 

95.266m 


De  tector 

y 


i 


-< 12.374m >] 


Fig.  4-21.  Ground  projection  of  field  of  view, 
based  on  the  cell  size  and  optics  as  shown  in  Fig.  4-22. 


Lense  System 
f.i.  = 1.7  cm 


Fig.  4-22.  Detector  cell  size  and  optical  system. 

As  a means  of  obtaining  some  feeling  for  the  length  of  time  that  a tar- 
get is  in  the  field  of  view  of  a given  cell,  and  the  fraction  of  the  total 
area  that  it  covers,  a computer  program  was  developed  to  plot  this  informa- 
tion. Figures  4-23  through  4-28  display  the  fractional  area,  6,  as  a func- 
tion of  time,  for  three  adjacent  cells.  Thus,  not  only  is  the  length  of 
time  the  target  is  in  any  one  cell  apparent,  but  the  time  displacement  from 
one  cell  to  the  next  is  evident.  This  information  is  relevant  when  it  is 


recalled  that  the  final  decision  operation  is  based  on  the  responses  of 
three  adjacent  cells. 
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Figure  ^-23 - Fractional  area  of  target  In  three  adjacent  cells. 


DELTA  FOR  MAN  VEL-5.0  ANG=180 


Figure  Fractional  area  of  target  In  three  adjacent  cells. 
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DELTA  FOR  DOG  VEL=5.0  ANG=135 
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Figure  4-26.  Fractional  area  of  target  in  three  adjacent  cells. 
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Figure  b-27.  Fractional  area  of  target  In  three  adjacent  cells. 
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DELTA  FOR  CAR  VEL=10  ANG=180 


Figure  l*-28.  Fractional  area  of  target  In  three  adjacent  cells. 
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Several  different  baclground  models  were  also  used  in  the  calculations. 


These  are  tabulated  in  Table  4-2. 


Table  4-2  Tabulation  of  Background  Models 


Background 

Temperature 

Emissivity 

Grass 

288 

.92 

Concrete 

316 

.95 

Soi  l 

288 

.39 

Composite 

294 

.72 

2.4.4.  Operation  at  Night — Se l ected  Results 

The  probability  of  detection  has  been  calculated  for  ni 
for  all  of  the  target  models  and  background  models  described 
Some  of  these  results  are  tabulated  in  Table  4-3.  The  first 
this  table  are  the  "worst  case"  models  used  to  determine  the 


ters. 


ght  operation 
in  Sec.  2.4.3. 
two  items  in 
system  parame- 
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Table  4-3.  Selected  Results  for  Night  Operation 


Target 

Background 

Angle  of 

Ent  ry 

Speed 

m/s 

PD 

Squi rrel 

Soil 

90° 

1 

5.6  x 10"8 

180° 

1 

6.9  x 10~10 

Man 

90°,  157.5°,  180° 

1 

1.0 

Crawl ing 

Grass 

157.5°,  180° 

.025 

1.0 

Black 

Car 

Any 

90°,  180° 

30 

1.0 

Red 

Car 

Any 

9C°,  180° 

30 

1.0 

o 

o 

o 

1 

1.0 

Dog 

Grass 

180° 

.025 

1.0 

J2.4.5  Operation  in  Day l ight — Se l ected  Resul ts 

The  probability  of  detection  has  been  calculated  for  daylight  operation 
(assuming  bright  sunlight)  for  all  of  the  target  models  and  background 
models.  Some  of  these  results  are  tabulated  in  Table  4-4. 
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Table  4-4.  Selected  Results  for  Day  Operation 


Angle  of 

Speed 

Background 

Entry 

m/s 

Squi rrel 

Soil 

90° 

180° 

1 

1 

2.54  x 10~2 

4.97  x 10~1A 

Man 

90°,  157.5°,  180° 

1 

1.0 

Crawl ing 

Grass 

157.5°,  180° 

.025 

1.0 

Black 

Car 

Any 

°o 

00 
* — 

o 

o 

o 

30 

1.0 

Red 

Car 

Any 

90°,  180° 

30 

1.0 

90° 

1 

1.0 

Dog 

Grass 

180° 

.025 

1.0 

Consideration  of  Atmospheric  Conditions 


The  atmospheric  attenuation  is  modeled  by  an  exponential  function 


t (R,X>=  e 


-a(X)R 


X = wavelength 


a(x)  = absorption  coefficient 


I 


R = range  of  observation 


It  has  been  observed  from  experirental  data  that  when  the  wavelength 
range  considered  is  in  the  3.4  um  to  4.2  um  region  then  the  atmospheric  at- 
tenuation is  nearly  constant  over  this  range  of  wavelength.  In  this  case, 
(4-1)  becomes 


T 


-n 


e 


1 


R 


(4-2) 


_ e-.2302aR 


where  = absorption  coefficient  in  (nepers/m) 

a = absorption  coefficient  in  (d  B/m) 

R = range  of  observation  in(m). 

■ 

Equation  (4-2)  was  used  in  the  calculations. 

The  values  of  the  absorption  coefficient  cover  a variety  of  atmospheric 

l 

conditions  ranging  from  a = .01  for  a clear  day,  to  a = .06  for  heavy  snow 
fall.  Calculations  have  been  made  for  absorption  coefficient  values  of 

a = .01,  .02,  .03,  .04,  .05,  .06,  .07,  (dB/m) 

The  effect  of  the  atmospheric  attenuation  is  to  reduce  the  signal-to- 
noise  ratio  of  the  target  models.  This  results  in  decreasing  the  probabili- 
ty of  detection. 

The  probability  of  detection  for  the  target  models  considered  has  been 
calculated  for  all  atmospheric  conditions.  However,  the  results  are  not 
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significantly  different  from  those  shown  in  Tables  4-3  and  4-4. 

2L.U.1_  Cons iderat  ion  of  Sun  Gl  int 

All  of  the  target  models  considered  thus  far  have  been  assumed  to  be 
diffuse  reflectors.  There  is  a possibility,  however,  that  motionless  ob- 
jects having  specular  reflecting  points  may  be  in  the  field  of  view  of  one 
or  more  cells.  Such  specular  reflecting  points  might  occur  on  parked  vehi- 
cles, aircraft,  or  other  objects  located  in  the  field  of  view  for  temporary 
storage.  Presumably,  objects  permanently  located  in  the  field  of  view  can 
be  painted  appropriately  to  eliminate  specular  reflecting  points. 

As  a result  of  the  sun's  motion,  it's  image  from  any  specular  surface 
may  appear  and  disappear  in  the  field  of  view  of  any  celt.  If  the  appeal — 
ance  or  disappearance  is  rapid  enough  the  resultant  change  in  illumination 
will  result  in  a decision  of  target  present.  Of  course,  once  the  image  has 
been  in  the  field  of  view  for  a period  of  time  longer  than  the  background 
averaging  time,  it  becomes  part  of  the  background  and  a decision  of  no  tar- 
get is  once  again  made. 

In  order  to  determine- the  seriousness  of  the  sun  glint  problem,  calcu- 
lations were  made  of  the  probability  of  detection  as  the  sun's  image  ap- 
peared in  a given  field  of  view.  For  the  purpose  of  these  calculations  the 
specularly  reflecting  surface  was  considered  to  be  a set  of  plane  facets 
such  that  the  sun's  image  can  enter  the  fields  of  view  of  three  adjacent 
cells  almost  simultaneously.  This  is  believed  to  be  a "worst  case"  assump- 
tion since  the  imagt  a curved  surface  would  be  more  diffuse  and,  hence, 
would  not  change  in  intensity  as  rapidly.  The  time  required  for  the  sun's 
image  to  cross  the  edge  of  the  surface  and  increase  to  maximum  intensity  is 
simply  the  time  required  for  the  sun  to  move  its  own  diameter.  Since  the 
sun  subtends  an  angle  of  0.533°  at  the  earth's  surface,  and  has  an  angular 
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velocity  of  .025°  per  minute,  the  time  in  question  is  about  128  seconds. 
The  background  averaging  time  for  the  proposed  system  is  40  seconds  so  it 
can  be  expected  that  significant  changes  in  image  intensity  can  take  place 
rapidly  enough  to  result  in  a decision  of  target  present. 

Calculations  of  the  probability  of  detection  have  been  made  assuming 
the  rate  of  change  implied  by  the  above  numbers.  In  all  cases  the  probabil- 
ity of  detection  reached  unity.  It  appears,  therefore,  that  sun  glint  may 
be  a potential  problem  in  the  proposed  system. 

One  way  of  alleviating  the  sun  glint  problem  is  to  reduce  the  back- 
ground averaging  time  so  it  becomes  part  of  the  background  before  the  target 
time  average  becomes  large  enough  to  result  in  a detection.  The  difficulty 
with  this  solution  is  that  it  would  also  reduce  the  probability  of  detection 
for  slowly  moving  targets  such  as  the  man  crawling  at  .025  meters  per 
second . 

^.4^.  Considerat  ion  of  Cl  oud  Motion 

Experimental  data  indicates  that  a cumulus  type  cloud  will  attenuate 
the  sun's  radiation  by  about  225  dB/km  of  cloud  thickness.  Thus,  a cloud 
drifting  between  the  sun  and  the  field  of  view  of  any  cell  will  drastically 
reduce  the  reflected  signal  from  the  background,  and  this  change  may  result 
in  a decision  of  target  present. 

In  calculating  the  probability  of  detection  due  to  cloud  motion,  a num- 
ber of  assumptions  were  made.  Specifically, 

1)  The  cloud  velocity  is  assumed  to  be  0.2  meters  per  second.  This 
corresponds  to  a slowly  moving  cloud  in  a breeze  of  about  .5  mph. 
Clouds  moving  faster  than  this  are  more  likely  to  be  detected. 
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2)  The  reflected  signal  from  the  background  is  assumed  to  be  attenuated 
linearly  with  time  until  it  reaches  zero  when  the  field  of  view  is  en- 


: 
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tirely  covered.  The  maximum  length  of  time  for  this  to  take  place  oc- 
curs when  the  cloud  is  moving  at  an  angle  of  9^  = C°  or  180"  and  is 
about  62  seconds.  The  change  would  be  more  rapid  for  other  angles  of 
approach  and,  hence,  detection  would  be  more  probable. 

Calculations  of  the  probability  of  detection  based  on  the  above  assump- 
tions yielded  a value  of  unity  even  when  the  emissivity  of  the  background  is 
0.92,  corresponding  to  the  minimum  reflected  signal.  Thus,  it  appears  that 
cloud  motion  may  also  be  a problem  with  the  proposed  system. 

Decreasing  the  background  averaging  time  is  not  likely  to  be  effective 
in  dealing  with  cloud  motion  because  under  conditions  of  higher  wind  veloci- 
ties the  clouds  will  be  moving  at  about  the  same  speed  as  many  of  the  tar- 
gets that  are  to  be  detected.  j 

A possible,  but  not  desirable,  solution  to  the  problem  of  cloud  motion 
is  to  simply  delay  declaring  any  decision  until  more  cells  have  been  affect- 
ed. From  this  information,  some  estimate  of  the  stze  of  the  affected  area 
can  be  made.  If  the  area  affected  is  larger  than  that  of  any  target  of  in- 
terest, the  disturbance  is  assumed  to  be  a cloud  and  a decision  of  no  target 
is  made.  However,  if  the  number  of  cells  affected  is  small  enough  to  indi- 
cate a target  of  normal  size,  a decision  of  target  present  is  made.  The  ob- 
vious problem  with  this  solution  is  the  large,  amount  of  delay  (on  the  order 
of  50  seconds)  required  to  announce  the  presence  of  any  target. 

2^. 9^.  Considerat  ion  of  Incandescent  Lamps 

The  presence  of  lighting  fixtures  in  the  field  view  of  any  cell  raises 
the  possibility  of  a false  alarm  resulting  from  the  light  being  turned  on  or 
off.  In  the  wavelength  region  of  interest  the  incandescent  lamp  is  the  one 
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that  is  most  likely  to  be  a problem  since  its  spectral  radiation  is  greater 
in  this  region  than  that  of  other  common  types  of  light  sources. 


I 
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For  purpose  of  calculation,  a 500  W lamp  operating  at  a temperature  of 
2960°  was  assumed.  It  was  further  assumed  that  80%  (i.e.,  400W)  of  the  to- 
tal input  power  is  radiated  with  a spectral  characteristic  corresponding  to 
black  body  emission  at  the  specified  temperature.  The  probability  of  detec- 
tion was  then  calculated  for  this  lamp  when  it  is  turned  on  or  off.  This 
probability  of  detection  in  each  cell  was  found  to  be  unity  in  all  cases. 
However,  if  there  is  a lamp  in  only  one  cell,  a target  will  not  be  declared 
because  the  adjacent  cells  are  not  affected.  If  there  are  lamps  in  two  or 
more  adjacent  cells,  being  turned  on  or  off  at  the  same  time,  then  a target 
will  be  dec l ared . 

There  are  two  obvious  ways  to  deal  with  this  problem: 

1)  Introduce  a shield  between  the  lamp  and  the  detector  that  is  opaque  in 
the  wavelength  range  of  the  detector.  It' may  or  may  not  be  opaque  in 
the  visible  range,  depending  upon  the  illumination  requirements.  With 
sufficient  mass,  or  sufficient  separation  from  the  lamp,  the  change  in 
temperature  of  this  shield  wili  be  slow  enough  to  avoid  detection. 

?)  Derive  signals  from  the  lamp  switching  system  to  instruct  the  detection 
system  to  ignore  the  observed  changes  until  they  have  been  absorbed 
into  the  background  average.  This  action  need  be  taken  only  for  the 
celts  whose  field  of  view  contain  the  lamp  so  that  all  other  cells  are 
in  full  operat ion. 

2.4.10.  Use  of  Optical  Fi l ters 

All  of  the  system  performance  reported  in  the  preceding  sections  as- 
sumed a wa'elength  range  in  the  optical  system  of  3.4  to  4.2  micrometers. 
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Some  sample  calculations  were  also  made  in  two  different  wavelength  ranges; 
one  range  being  from  3.2  to  3.54  micrometers  and  the  other  from  3.75  to  4.2 
micrometers.  In  both  cases  it  is  assumed  that  cold  filters  are  used  so  that 
their  emission  in  the  stop  bands  can  be  ignored. 

In  both  wavelength  regions  the  probability  of  detection  is  essentially 
the  same  as  it  is  for  the  original  band  for  the  same  class  of  targets  and 
background.  The  major  difference  is  that  saturation  occurs  in  the  3.2  to 
3.54  mm  range  for  stare  times  of  0.1  second  because  of  the  increased  detec- 
tor response  in  this  wavelength  interval.  With  the  use  of  stare-time  con- 
trol this  is  not  a problem,  but  there  does  not  appear  to  be  any  advantage  to 
using  the  smaller  wavelength  interval  since  system  noise  is  not  a limiting 
factor. 
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2.5.  SYSTEM  IMPLEMENTATION 


2.5^.  1 . Microprocessor  Design  for  £ Single  Array 

In  designing  a microprocessor  system,  one  major  concern  is  the  speed  of 
execution  of  the  computation.  In  the  proposed  system  about  0.1  second  is 
available  between  the  collection  of  data  sets  for  computations.  Some  mi- 
croprocessors, such  as  the  popular  Intel  8080  and  Motorola  6800,  are  rela- 
tively slow,  with  clock  cycles  of  1 microsecond.  Also, one  instruction  may 
require  more  than  one  clock  cycle  to  execute.  This  application  requires 
much  data  manipulation  for  vectors  of  length  256,  and  these  popular  M0S  pro- 
cessors are  too  slow  for  the  proposed  system.  For  example,  if  an  instruc- 
tion must  be  applied  to  all  256  data  words  of  the  vector,  and  that  instruc- 
tion requires  10  microseconds  to  execute,  then  the  one  operation  would  re- 
quire 2.56  milliseconds  to  complete.  This,  of  course,  does  not  consider  all 
the  overhead  involved  in  order  to  get  ready  to  execute  the  instruction,  such 
as  bringing  the  data  from  memory  into  the  processor. 

Some  microprocessors  are  available  with  clock  cycle  times  of  100  ns. 
They  are  bipolar  TTL  and  each  chip  has  fewer  transistors  than  the  M0S  pro- 
cessors described  above.  Thus,  they  are  functionally  less  complex,  being 
able  to  process  only  2 or  4 bits  of  information  per  chip  as  compared  to  8 
bits  for  the  8080  and  the  6800.  Also,  these  bipolar  chips  are  "micropro- 
grammed". That  is,  they  have  no  instructions  defined  for  the  processor. 
External  circuitry  must  define  the  instructions  and  supply  the  control  sig- 
nals to  the  processor.  It  is  desirable  to  have  more  than  2 or  4 bits  in  a 
computer  and,  fortunately,  these  microprogrammabl e processors  can  be  com- 
bined to  provide  longer  word  lengths.  The  Texas  Instrumens  SN54/74S481  is  a 
4 bit  processor,  and  two  chips  can  be  combined  to  make  one  8 bit  processor. 
These  "bit-slice"  processors  require  more  circuitry  to  control  them,  and 
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more  effort  to  design  a workable  system. 

For  the  present  application  it  is  desirable  to  have  a processor  with 
the  speed  of  the  bipolar  devices  without  spending  the  time  to  design  all  the 
instructions.  Such  an  alternative  is  available  in  the  Signetics 
300KT8080SK,  an  8080A  emulator.  This  processor  fits  on  one  printed  circuit 
board  and  is  provided  with  all  the  same  instructions  as  an  Intel  8080A  pro- 
cessor. In  addition,  12  instructions  are  allowed  to  be  microprogrammed  by 
the  system  designer.  For  instance,  a square  root  function  might  be  needed 
for  a certain  catcultion.  Rather  than  make  the  programmer  write  the  code 
for  the  algorithm,  an  instruction  to  perform  the  algorithm  can  be  written 
into  the  processor.  The  emulator  has  a cycle  time  of  100  ns,  and  it  is 
claimed  that  the  speed  of  operation  is  from  2 to  9.2  times  faster  than  the 
8080A  that  it  emulates. 

This  8080A  emulator  is  a good  preliminary  choice  for  the  IRCCD  Intru- 
sion Detection  System.  In  the  following  pages,  a program  is  designed  and 
the  memory  and  time  needs  determined.  It  is  shown  that  the  processor 
described  above  can  meet  the  time  constraints. 

The  flow  chart  of  Fig.  5-1  shows  the  general  scheme  for  computing  an 
alarm  decision.  Flow  charts  in  Figs.  5-2  through  5-5  elaborate  on  each 
piece  of  the  scheme. 

The  program  must  first  be  initialized,  as  in  Fig.  5-2.  Certain  con- 
straints are  needed  for  the  main  routine  to  operate.  Also,  a beginning 
value  of  m(t')  is  needed.  Figure  5-2  suggests  that  N data  sets  and  their 
sums  be  taken  in  order  to  make  this  computation.  An  alternative  method  is 
to  provide  initial  values  for  m(l')  and  the  cummulative  sums  used  to  compute 
it.  Once  the  initial  data  are  ready,  NOW  is  set  to  point  to  the  address  of 
the  newest  data,  and  the  main  program  begins. 


Figure  5~1.  General  flow<:  tation  of  IR  detection 

system. 


Figure  5-3  shows  the  steps  for  computing  mU').  First,  check  to  see  if 
m(t')  should  be  computed  now  by  checking  if  PCK  = P.  If  it  is,  continue 


with  the  calculation.  Otherwise,  load  new  data,  and  continue  the  program. 
The  background  time  average  is  defined  as 


m(i.')  = (1/N)  23  xCPU'-nn 

n = 1 


(5-1) 


If  mCt')  were  computed  directly  for  each  sensor  i,  then  each  time  that  it 
was  computed  (N-1)  additions  would  be  required.  Clearly,  for  large  N,  this 
requires  an  unacceptable  amount  of  time. 

An  alternative  is  to  store  cummulative  sums  for  m(n'),  where 


MSUM(i)  = 23  xCPU'-m,i)] 

n = 2 


(5-2) 


If  x(j,',i)  is  an  eight  bit  number,  then  MSUM  wilt  be  at  most  a 16  bit  numer 
if  N = 16.  For  each  new  data  set,  m(s.')  can  be  computed  by  adding  the  new 
data  and  subtracting  the  oldest  data  from  the  cummulative  sum.  So, 

MSUM(i)  = MSUM(i)  + new  data  / 

MU')  = (1/N)  MSUM(i) 

MSUM(i)  = MSUM(i)  - oldest  data.  f 

To  find  the  oldest  data,  a list  ordered  by  age  is  kept  of  the  addresses 
of  all  the  data  used  in  computing  m(t').  As  new  data  arrived,  the  start-, 
address  of  the  data  will  be  added  to  the  end  of  the  queue,  and  the 
data  address  will  be  deleted  from  the  front  of  the  queue.  An  explana' 
queues  and  their  uses  is  found  in  "The  Art  of  Computer  Programming  , 
by  D.  E.  Knuth,  Addison-Westey,  1973,  pages  240-242. 
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F?gi;re  5-3.  Flowchart  for  computation  of  m(jT),  the  background  time 
average. 


The  method  above  requires  a small  amount  more  memory  space  than  the 
method  that  computes  the  entire  sum  each  time  new  data  arrives.  However,  it 
requires  significantly  less  time  to  compute. 

To  compute  ^m(t')  the  fastest  method  is  table  look-up.  If  mU')  is  8 
bits  long,  then  all  'a»jes  of  *7  are  stored  in  a table  in  256  memory 
locations.  If  the  table  begins  at  location  A,  then  is  located  at 
location  A + mU'). 

Next,  a new  data  set  is  loaded  and  its  spatial  average  computed,  as  in 
Fig.  5-4.  It  is  recommended  that  some  sort  of  direct  memory  access  be 
available  for  the  incoming  data  to  be  loaded  into  memory.  A program  con- 
trolled data  input  would  not  be  fast  enough  to  load  one  datum  each  4 mi- 
croseconds. D should  be  a power  of  2,  for  then  division  may  be  accomplished 
by  a series  of  right  shifts,  one  for  each  factor  of  2.  If  the  spatial  aver- 
age is  too  high,  t is  shortened;  if  it  is  too  low,  t is  lengthened. 

Finally,  the  threshold  and  alarm  computations  are  made,  as  in  Fig.  5-5. 
First,  wCD  is  computed  for  all  256  sensors.  The  alarm  computations  can  be 
computed  in  two  steps.  First,  v(i),  a vector  of  length  254,  is  computed, 
where  each  v(l,i)  is  a 3 bit  word  and  reflects  a comparison  of  b(*,i)  with 
w(l,i),  w(l,i+1),  and  w(l,i+2).  Second,  v(D  is  scanned  n/3  words  at  a time 
to  see  if  t or  more  l 1 s are  present  in  the  group.  But,  computing  v(t)  is 
unnecessary  and  is  bypassed.  Fig.  5-5  shows  a computation  for  n * 12  and  t 
= 8.  For  each  set  of  comparisons  between  b(t,i)  and  the  three  consecutive 
channels  of  w(l),  a count  of  the  number  of  I ' s is  kept.  When  n/3*4  counts 
have  been  made,  the  sum  of  the  4 counts  is  checked  to  see  if  it  is  greater 
than  8.  If  it  is,  an  alarm  decision  is  made.  If  the  count  is  less  than  8, 
the  search  for  an  alarm  decision  is  resumed  by  eliminating  the  oldest  count, 
c(0),  and  computing  the  next  count,  c<3>. 
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This  program  can  be  divided  into  four  main  parts:  the  computation  of 
m(n'),  the  computation  of  the  spatial  average,  the  computation  of  wU),  and 
the  threshold  and  alarm  decision.  For  each  of  these  four  pieces,  an  esti- 
mate of  the  running  time  is  based  on  the  execution  time  of  the  instructions 
for  the  8080A  emulator  (Signetics  Corporation,  "Technology  Leadership  Bipo- 
lar Microprocessor",  Feb.  1977,  pages  42-43).  These  numbers  are  intended  as 
ballpark  figures  and  are  not  absolute  ones.  The  actual  running  time  will 
depend  on  the  cunning  of  the  programmer. 

The  computation  time  estimate  is  as  follows: 

Compute  m ( t. ' ) : 

load/store  1592(0.9)  = 1432.8  microseconds 

add/subtract  512(1.05)  = 537.6 

shifts  (for  division  for  N=32) 

1280(0.6)  = 668.0 
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Compute  ) : 

double  precision  load 
256(1 .8) 

add  256(.6) 

double  precision  load 
256(1.8) 
multiply  256(3.9) 


Compute  Spatial 

Average 

load 

18(0.9) 

add 

17(0.6) 

shift 

5(0.6) 

compare 

2(0.6) 

jump 

2(0.9) 

output  t 

1(1.2) 

= 460.8 

= 153.6 

= 460.8 

= 998.4 


4712.0  microseconds 

* 16.2 

* 10.2 

* 3.0 

« 1.2 

* 1.8 

« 1.2 


33.6  microseconds 

93 


r 

pi 

Compute  w(i,i):  ; 


toad 

256(0.9)  = 

230.4 

subtract 

256(1.05)  = 

268.8 

compare 

256(0.6)  = 

153.0 

jump 

256(0.9)  = 

230.4 

compl ement 

256(0.3)  = 

76.8 

increment 

256 (.45)  = 

115.2 

\ 

store 

256(0.9)  = 

230.4 

J 

* 

J 

1305.0  microseconds 

w 

> 

Compute  Threshold  and  Alarm 

• 

• 

t 

compare 

1266(0.6)  = 

759.6 

l 

' 

jump 

1266(0.9)  = 

1139.4 

add 

754(0.6)  = 

452.4 

J 

toad 

1016(0.9)  = 

914.4 

* j 

3265.8  microseconds 


Total  Computation  Time 


I 

t 


Total  9.316  ms 

The  Random  Access  Memory  (RAM)  needed  for  data  storage  may  be  estimated 
as  follows: 

data  sets: 

for  N = 32:  256  x 8 bits  = 8 K words 

m(l):  256  x 8 bits  = 1/4  K words 

MSUM:  256  x 16  bits  = 1/2  K words 

w(t):  256  x 8 bits  = 1/4  K words 

b(D:  256  x 8 bits  = 1/4  K words 

9 1/4  K words 

These  are  the  five  largest  users  of  storage.  Let  us  estimate  that  12 
K,  8 bit  words  of  RAM  be  used.  This  estimate  is  based  partially  on  the  fact 
that  often  RAM  may  be  bought  in  multiples  of  4K  words. 
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2.5.2.  Est imated  Cost  for  £ Single  Array 

The  system  cost  and  power  need  may  be  estimated  as  follows: 

Estimated  Price  Estimated  Power 


8080A  emulator: 

Signet ics  300KT8080K  $ 300  4 amps 

4K  ROM  for  program  storage 

4 82S184's  160  1 amp 

12K  RAM  for  data  storage 

96  S4015-3's  600  10  amps 


$1060  15  amps  at  5 volts 

In  the  future,  new  products  may  make  these  choices  of  ROM  and  RAM  ob- 
solete. If  the  system  is  built,  these  choices  should  be  reevaluated. 

A typical  power  supply  choice  might  be  DATEL's  MPS-5/18,  a 5 volte,  18 
amp  supply  with  a cost  of  about  $100. 

Packaging  costs  must  be  added  to  this  parts  estimate.  These  should  be 
determined  by  the  environmental  needs  when  the  system  is  built. 

2.5.3.  Time-Sharing  for  Mul tipi e Arrays 

It  may  be  possible  to  use  one  processor  to  control  more  than  one  sensor 
array.  However,  the  number  of  arrays  that  may  be  controlled  will  be  limited 
more  by  the  amount  of  memory  needed  than  by  the  time  available  for  computa- 
tions. The  8080A  emulator  described  above  may  directly  access  up  to  64K 


words  of  memory.  If  each  array  requires  10K  words  for  data  storage,  then  at 
most  64/10  - 6 arrays  may  be  serviced  by  one  processor  due  to  the  direct 
memory  accessing  limitations  of  the  processor.  If  our  ballpark  figure  for 
the  computation  time  of  our  program  is  reasonable,  then  6 arrays  could  be 
serviced  in  the  given  0.1  seconds. 

On  the  basis  of  the  cost  and  power  estimates  given  above,  the  cost  for 
6 arrays  would  be  on  the  order  of  $4460  (not  including  power  supplies  and 
packaging)  and  the  power  requirements  would  be  61  amperes  at  5 volts.  Thus, 
the  cost  per  array  will  be  on  the  order  of  $750. 


3.1.  CONCLUSIONS 


Based  on  the  analysis  presented  in  this  report,  and  upon  the  computa- 
tions that  have  been  performed,  there  are  a number  of  conclusions  that  ap- 
pear to  be  justified.  Most  of  these  have  been  stated  in  earlier  sections  of 
this  report,  but  they  are  restated  here  in  a more  concise  form  as  a matter 
of  collecting  them  together  in  a single  place. 

1)  There  is  strong  analytical  and  computational  evidence  that  the  IRCCD 
array  can  be  successfully  used  for  intrusion  detection. 

2)  By  employing  a decision  rule  that  utilizes  the  outcomes  of  multiple  ob- 

servation (in  both  space  and  time)  the  detection  characteristic  can  be 
made  sufficiently  steep  to  minimize  false  alarms  due  to  small  targets 
while  maximizing  the  probability  of  detection  for  targets  of  interest. 

3)  False  alarms  due  to  system  noise  appear  to  be  negligibly  small  when  the 

computed  thresholds  are  used. 

4)  The  most  serious  problem  appears  to  be  false  alarms  due  to  sun  glint 

and  cloud  motion. 

5)  It  appears  quite  possible  to  implement  the  electronics  required  for 

performing  the  decision  function  with  a single  microprocessor  and  some 
associated  memory.  The  cost  of  this  elections,  at  present  prices,  is 
probably  under  $1200,  not  including  packaging. 

6)  It  also  appears  possible  to  use  a single  microprocessor,  with  addition- 
al memory,  to  perform  the  decision  function  for  as  many  as  six  IRCCD 
arrays.  The  cost  of  the  electronics,  per  array,  is  probably  under 
$900,  not  including  packaging. 
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3.2.  RECOMMENDATIONS 

The  principal  recommendation  is  that  the  next  step  in  the  development 
of  a working  system  be  taken.  This  recommei  dat ion  is  based  on  the  analyti- 
cal and  computational  results  that  have  been  obtained  in  this  study.  It  is 
believed  that  these  results  are  sufficiently  favorable  to  indicate  a high 
probability  that  a successful  system  can  be  developed,  with  the  system  pro- 
posed in  this  report  serving  as  a reasonable  point  of  departure  in  this 
development.  The  next  logical  step  is  to  construct  a breadboard  model  of 
system  ana  test  it  under  as  realistic  conditions  as  possible. 

If  a decision  is  made  to  develop  the  proposed  system  through  the  bread- 
board stage,  there  are  several  analytical  and  computational  tasks  tnat  need 
to  accompany  the  more  detailed  system  design.  The  following  recommended 
tasks  are  approximately  in  the  order  in  which  they  should  be  done,  although 
of  course  some  of  them  can  be  carried  out  in  parallel. 

1)  Review  the  available  microprocessors  and  memories  with  respect 
to  the  system  requirements,  select  tne  components  needed  and 
order . 

2)  Obtain  more  accurate  information  regarding  the  IRCCD's  avail- 
able at  that  time  and  re-evaluate  the  system  analysis  with 
respect  to  selecting  firm  values  of  stare  time,  background 
averaging  time,  threshold  levels,  and  observation  multiplici- 
ty. 

3)  Write  the  programs  necessary  to  accomplish  the  desired  compu- 
tations, after  checking  the  flowcharts  to  establish  that  they 
are  correct  and  complete. 
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4)  Prepare  a test  program  designed  to  establish  the  performance 


of  the  system  relative  to  the  design  objective. 
5)  Construct,  de-bug  and  test. 
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APPENDIX  4 

OPTIMUM  36T5CTI0N  SYS'E* 

The  analysis  presented  here,  in  abbreviated  fan,  follows  the  genera, 
method  outlined  in  C2,  Chao.  2-  It  differs  from  that  analysis  only  in  that 
the  noise  level  here  is  a function  of  signal  amplitude  rather  tnar  be^ng 
constant . 

It  is  assumed  that  M observations  from  each  of  K detector  cells  are 
available.  Thse  observations  form  the  elements  of  a vector. 

r s (r^ ,r . • .,r^)  , n = WK  (A-1) 

Target  detection  is  based  on  a likelihood  ratio  test  defined  by 


t(r) 


pCrjH.j) 

pTTTTT) 


> 

< 


x 


(A-2) 


where  HQ  is  a hypothesis  that  no  target  is  present  in  any  of  the  HK  observa- 
tions, while  H1  is  the  hypothesis  that  a target  is  present  in  all  observa- 
tions. The  decision  threshold,  x,  is  determined  by  the  desired  probability 
of  false  alarm  and  evaluated  by  the  MK-fold  integral 

P.  * / p(r|H  ) dr_  (A-3) 

f Hr)  > x 

The  probability  of  detection  yielded  by  this  threshold  is 

Pn  = / p(r |H«)  dr  (A-4) 

Hr)  > x “ 

The  vector  jr  is  assumed  to  be  multivariate  Gaussian  on  the  basis  of  the 
Gaussian  approximation  to  the  Poisson  distribution  for  the  output  of  each 
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cell.  Thus,  the  conditional  density  function  p(rjHQ)  is 

p(r|H  ) = (2ir)"n/2|A  |~1/2  exp  [-1/2(r-m  )A_1  (r-m  )T:  (A-5) 

where 

% =(mo1'  mo2'"-'mon) 
moi  = ECrilV 

ao1n 
aonn_ 


ao11  a012 

• • • 

a * • • 

onl 


01] 


= a 


oji 


ECC(r.-m  .)|H  ][r.-m  .)|H  D> 
1 01  1 0 ] 0]  1 o 


and  the  superscript  T denotes  a transpose.  It  is  assumed  here  that  (r.-m  .) 

1 01 

and  (rj-moj)  are  uncorrelated  for  i t j since  they  represent  outputs  from 
different  cells  or  at  different  times.  Thus, 


i * j 


In  a similar  way  the  conditional  densit,  function  p(rjH.j)  is 

P(£.l H1 ) = (2w)"n/2  |A1f1/2  exp  C-1/2(£-m1)A^1(£-m1)T]  (A-6) 

where 
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— 1 = (*11,  m12,  m1n) 


"H  = ECri'H1] 


a1 11  a1 12 


'lij  = a^.  *E«(r.-m1i)|H1K(rj-n1j)+H1]> 


and , as  before. 


am  * Or  i * j 


(A-7) 


i = J 


Using  (A-5)  and  (A-6),  the  likelihood  ratio  becomes 


- jr^rny  - <■ 


1/2C(£-£^ ) 1 (^-n^  )T 


C A— 8) 


- (r-m  > A”1 (r-m) T]> 
— 0-0 0 


It  is  more  convenient  to  consider  the  natural  logarithm  of  (A-8)  and  write 


it  as 


inCKr,)]  - = 1/2[(r_-mQ)  A^Cr-n^)7  - (r-m,, ) A^1  ir-m,)1]  (A-9) 


Hence,  a sufficient  statistic  can  be  defined  as 


i 


L(r)  = 2 


tn[l(r)J  - l 


(A-10) 


from  which  the  likelihood  ratio  test  becomes 


H 


1 


L(r)  = ( r -m  ) A 1 (r-m  )T  - (r-m,)  Ai^  (r-m.,)T  > n 
— o — o — o i—i  — — i < 

H 


(A-11) 


and 


Pc  = / p(r |H  )dr  (A-1 2) 

f L(r»n  ~ 0 ~ 


and 


pn  = / p(£.|H1)d£  (A-1 3) 

L(^»n  1 

Because  the  covariance  matrices  Ao  and  A^  are  different  under  the  two 
hypotheses,  it  is  not  possible  to  reduce  the  decision  boundary  defined  by 
(A-11)  to  a plane  in  n-dimensiona(  space.  Thus,  further  analytical  evalua- 
tion of  (A-12)  and  (A-13)  is  not  possible  and  their  evaluation  by  computer 
is  extremely  cumbersome  when  n is  greater  than  2.  However,  in  order  to  com- 
pare the  theoretical  optimum  results  with  those  obtained  by  the  proposed 
system,  it  is  desirable  to  carry  out  such  an  evaluation  for  at  least  the 
simplest  case. 

The  simplest  case  corresponds  to  M=K*na1.  In  this  case 
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T' 


p 


“ r 
i»  * nr 

— o o 
■l  * *1 

A * of  , A"1  * i/02 
— O 0 —0  o 

A,  * of  , A,'1  * 1/of 

Thus, 


Kr)  * <r-mo>  (-iy)  (r-mo)  - (r-m^ ) (r-»^ ) 

°o  °1 


(A-H) 


In  the  present  case,  the  variances  are  proportional  to  the  mean  values 
because  of  the  properties  of  the  underlying  Poisson  distribution.  For  this 
situation,  the  coefficient  of  r in  (A-14)  vanishes  and  the  likelihood  ratio 
test  becomes 


2 > Vi  r 

r < m -m  n 

Hi  1 ° 


(m  - m. ) 3 = n 1 
o 


In  either  case 


Pp  = f ( 77oo)  **  exp[-|V-m^j /2o^3dr 


( A— 1 7 


r2>n* 


Pjj  = f ( 7To^)  ^ expC-jj"-m J2/2o2]dr 


(A-18 


r2>n 1 


Hence,  n'  is  determined  by  selecting  a desired  value  of  Pp,  and  from  this 


P0  can  be  calculated. 
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